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^This  final  report  describes  all  the  projects  undertaken  during  the  entire 
period  of  the  above  contract,  the  scope  of  which  was  to  study  coronal  holes  at 
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previously  published. 
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in  the  Stanford  9.1  cm  synoptic  charts  during  the  same  period,  when  clear  X-ray 
data  of  coronal  holes  were  available  from  Skylab. 

Chapter  2  summarizes  our  efforts  to  identify  Individual  high  latitude /Corona] 
holes  with  high  speed  solar  wind  streams  far  above  or  below  the  ecliptic,/  The 
coronal  holes  were  Identified  from  the  Kltt  Peak  10830  X  synoptic  maps,  while  the 
high  speed  solar  wind  streams  were  Identified  from  the  Interplanetary  scintilla¬ 
tion  measurements  made  at  UCSD  and  at  the  Toyokawa  Observatory  in  Japan. 

Chapter  3  describes  our  1977  coronal  hole  observations  with  the  Haystack 
(1.94  and  3.81  cm)  and  Arecibo  (11.5  and  21.4  cm)  radio  telescopes./  The  Arecibo 
observations  were  of  special  importance  because  they  demonstrated  tfiat  in  the 
radio  domain  coronal  holes  become  most  clearly  visible  in  the  10-10{0  cm  range, 
provided  the  radio  telescope  has  sufficient  angular  resolution.  These  observa¬ 
tions  produced  a  radio  spectrum  of  a  coronal  hole. 

Chapter  4  describes  the  international  effort  we  organized  in  the  summer  of 
1978  to  obtain  solar  radio  maps  over  a  wide  range  of  wavelengths.  These  radio 
maps  were  made  on  the  same  days  at  a  number  of  radio  observatories  around  the 
world.  ']  During  the  first  observing  period,  July  6-10,  1978  there  were  no  coronal 
holes  present.  During  the  second  period,  July  15-18,  1978,  there  were  two  siz¬ 
able  coronal  holes,  one  in  the  northern  hemisphere  and  the  other  in  the  southern 
hemisphere.  Unfortunately  there  were  also  several  active  regions,  which  made  it 
difficult  to  establish  the  quiet  sun  level  and  Introduced  some  uncertainty  in  the 
radio  spectra  obtained  for  these  two  coronal  holes. 

•  Chapter  5  describes  the  use  of  the  observed  radio  spectra  to  produce  models 
of  the  solar  atmosphere  inside  coronal  holes.  A  ray  tracing  program  is  used  in 
conjunction  with  a  standard  model  for  the  quiet  solar  atmosphere  to  compute 
brightness  temperatures  at  different  radio  wavelengths,  initially  for  the  stan¬ 
dard  model  and  then  for  a  series  of  modifications  affecting  the  temperature, 
electron  density  and  neutral  particle  density  profiles  of  the  standard  model  in 
an  effort  to  find  the  profiles  that  would  reproduce  the  radio  spectra  obtained 
from  our  observations. 

The  results  of  this  analysis  indicate  that  inside  coronal  holes  the  follow¬ 
ing  changes  occur:  The  chromosphere  expands  in  height  by  about  a  factor  of  2, 
while  the  temperatures  and  electron  and  neutral  particle  densities  decrease  to 
about  80%  of  their  quiet  sun  values.  The  critical  transition  zone  expands  by  a 
factor  of  about  5-6,  while  the  temperatures  decrease* to  about  50%  and  the  elec¬ 
tron  densities  to  about  20%  of  their  quiet  sun  values./  Finally,  the  corona 
shows  no  substantial  expansion,  but  the  temperatures  and  electron  densities 
decrease  to  about  80%  of  their  quiet  sun  values. 


Appendix  1  lists  the  model  for  the  quiet  solar  atmosphere  used  in  this 
study.  Appendix  2  gives  a  complete  listing  of  the  computer  program  used  in  the 
computation  of  the  radio  spectra  for  modified  atmospheric  profiles,  as  well  as 
a  sample  input  and  output. 
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INTRODUCTION 


I .  Scientific  Objectives 

The  purpose  of  this  contract  was  to  study  coronal  holes 
at  radio  wavelengths  and  to  use  these  radio  observations  to: 
1)  Identify  coronal  holes  in  solar  radio  maps  at  different 
wavelengths,  2)  model  the  structure  of  a  coronal  hole,  3) 
compare  different  coronal  holes  in  order  to  see  If  there 
were  observable  differences  at  radio  wave  lengths ,  and  4) 
correlate  the  coronal  hole  observations  with  high  speed 
solar  wind  streams. 

The  structure  of  coronal  holes  is  still  not  well 
understood.  Efforts  to  model  coronal  holes  have  been  made  by 
several  authors  (Withbroe  and  Wang,  1972;  Mariska,  1973; 
Mariska  and  Withbroe,  1978)  using  EUV  data.  In  addition, 
Chluderi-Drago  (  Ch  iud  e  r  i-I)  rag  o  ,  et  al.,  1  977  ;  Ch  iuder  i-Drago 
and  Poletto,  1  977  ;  Ch iude r  i- D r ago  ,  1980)  has  made  several 
attempts  to  model  coronal  holes  based  on  radio  observations. 
Finally,  Dulk  and  his  co-workers  (Dulk,  et  al.,  1977)  have 
attempted  to  match  the  models  deduced  from  EUV  and  radio 
data.  Unfortunately,  the  models  derived  from  EUV  data  are  in 
substantial  disagreement  with  those  derived  from  radio  data. 
It  is  Important  to  resolve  these  differences  and  much  of  our 
work  has  been  aimed  toward  this  problem. 

A  clear  understanding  of  the  transition  zone  is  of  great 
importance  to  our  understanding  of  the  coronal  hole 
phenomenon.  As  Withbroe  (1978)  has  pointed  out,  "The 


6 


transport  of  energy  by  thermal  conduction  into  the 
transition  zone  is  a  primary  coronal  cooling  mechanism." 
Thus,  the  temperature  and  particle  density  gradients  in  the 
transition  zone  are  of  critical  importance  in  understanding 
the  energy  balance  of  a  coronal  hole. 

The  link  between  large  equatorial  coronal  holes  and  high 
speed  solar  wind  streams  observed  at  the  earth  has  been  well 
established  (Mewpert  and  P  i  z  z  o  ,  1976;  Molte  et  a  1  . ,  1976, 
1978;  ilroussard,  et  al.  ,  1  9  78).  Smaller  coronal  holes  and 
coronal  holes  at  high  solar  latitudes  have  also  been  related 
to  high  speed  solar  wind  streams  (Coles  and  Rickett,  1976; 
Kickett,  et  al.  1976;  Sine  and  Rickett,  1973;  Baker  and 
Papagiannls,  1980),  although  the  correlation  is  less  clear 
for  these  cases.  A  qualitative  understanding  of  the 
relationship  between  a  coronal  hole  and  high  speed  wind 
streams  emanating  from  the  coronal  hole  has  been  provided  by 
Konp  and  Orrall  (1977),  but  a  quantitative  understanding 
requires  a  detailed  knowledge  of  the  energy  balance  in  the 
coronal  hole.  Thus  we  are  again  brought  back  to  the 
conclusion  that  consistent  models  of  coronal  holes  are 
required  to  complete  our  understanding  of  the  solar  wind 
streams. 

This  final  report  provides  a  complete  summary  of  our 
work  on  coronal  holes,  using  radio  observations  as  well  as 
interplanetary  scintillation  and  coronal  hole  maps  made  from 
RIJV  data,  X-Ray  data,  and  the  Helium  10330  ^ 
spectroheliogrnms  made  at  Kitt  Peak.  Chapter  1  describes  the 
identification  of  coronal  holes  on  the  synoptic  solar  radio 
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maps  made  for  other  purposes  at  Stanford  during  the  Skylab 
period  In  1973,  as  well  as  the  determination  of  a  coronal 
hole  radio  spectrum  from  different  observations  made  during 
the  same  period.  Chapter  2  describes  the  work  we  have  done 
on  the  correlation  between  high  latitude  coronal  holes  and 
high  speed  solar  wind  streams  observed  far  above  or  below 
the  ecliptic  by  Interplanetary  scintillation  measurements. 

Chapter  3  is  devoted  to  our  observations  made  at  the 
Haystack  Radio  Observatory  and  at  the  Arecibo  Radio 
Observatory  during  1977.  The  work  described  in  Chapters  1  - 
3  has  been  described  in  detail  in  our  three  previous 
Scientific  Reports  ( AFGL-TR-77-0292  ,  AFGL-TR-73-023 A , 
AFGL-TR- 3 0-0123) . 

Chapter  4  is  devoted  to  the  presentation  and  analysis  of 
the  solar  radio  maps  made  during  the  month  of  July,  1978. 
These  maps  were  made  as  part  of  an  international 
co-operative  effort  organized  by  us  to  provide  solar  maps 
over  a  wide  range  of  radio  wavelengths.  Finally,  Chapter  5 
describes  the  modeling  process  we  have  used  in  analyzing  the 
radio  spectra  we  have  obtained.  This  modeling  process 
provides  us  with  a  vertical  profile  of  temperature  and 
electron  density  inside  a  coronal  hole. 

In  addition,  there  are  two  appendices.  Appendix  I  gives 
the  electron  density,  neutral  particle  density,  and 
temperature  profiles  of  the  photosphere,  chromosphere, 
transition  zone  and  corona  that  we  have  used  as  a  standard 
model  of  the  quiet  solar  atmosphere.  Appendix  IT  gives  a 
listing  of  the  computer  program  that  was  used  to  modify  the 
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st  and  aril  solar  atmosphere  model  and  then  trace  ray  paths 
through  the  atmosphere  to  determine  radio  brightness 
temperatures  at  different  wavelengths.  In  addition  to  the 
program  itself,  we  have  included  a  sample  input  and  output 
of  this  p  rog  r am . 

2 .  Personnel 

The  period  covered  by  this  contract  extended  from  July 
1,  1977  to  September  30,  1980.  The  principal  investigator 

was  'Prof.  Michael  0.  Papagiannis.  From  July,  1977  to 
September,  1978,  Or.  Fred  Wefer  worked  on  the  project  as  a 
post-doctoral  research  associate,  while  from  January  1, 
1979  to  the  end  of  the  contract  period,  this  position  was 
held  by  Or.  Kile  Baker.  At  various  times  during  the  period 
covered  by  the  contract,  three  students  (Terry  Varner, 
Michael  Van  Steenberg,  and  Michele  de  la  Pena)  worked  on 
various  tasks  of  the  project. 

3  .  Puhl icat ions 

During  the  period  covered  by  the  contract,  two  papers, 
four  scientific  reports  including  this  final  report,  and 
five  abstracts  have  been  published  pertaining  to  the  work  of 
the  contract.  These  publications  are  listed  below: 

Papers 

"Coronal  Holes  at  11.5  and  21  cm  observed  with  the 
Areclbo  Radio  Telescope",  M.D.  Papagiannis  and 
F.L.  liefer;  Mature,  2_7_3>  5  20  ,  1  978. 


"Synoptic  Charts  of  Solar  9.1  cm  and  Coronal  Hole 
Data",  F.L.  Wefer  and  M.D.  Pa  pan,  Iannis;  Solar 
Physics,  67_,  1  3  ,  1  980. 

Scientific  Reports 

"The  Radio  Spectrum  of  Coronal  Hole  1",  F.L.  liefer  and 
M.D.  Papagiannis;  AKGL-TR-  7  7-029  2  ,  Boston  Unlv.  , 
Ser.  II,  No.  65,  1977. 

"Synoptic  Charts  of  Coronal  Holes  and  9.1  cm  Solar  Data 
from  May  Through  August  1973",  F.L.  Wefer,  M.D. 
Papagiannis,  M.K.  Van  Steenberg,  T.M.  Varner; 
AFGL-T  I? -78-0234,  Boston  llniv.,  Ser.  II,  No.  6  3, 
1  9  79  . 

"Correlation  of  High  Latitude  Coronal  Holes  with  Solar 
Wind  Streams  Hi? h  Above  or  Below  the  Ecliptic", 
K.B.  Baker  ana  M.D.  Papagiannis;  AFGL-TR-80-0 1 2 3 , 
Boston  Unlv. ,  Ser.  II,  No.  71,  1980. 

"Radio  Studies  o  Coronal  Holes",  M.l).  Papagiannis  and 
K.B.  Baker;  A  FG  L- TR- 3  I  -  00  5  3  ,  Boston  llniv.,  Ser. 
II,  No.  7  2  ,  1  98  1  . 

Abs  tracts 

"Coronal  Hole  Observations  at  1420  and  2600  MHz  with 
the  Areclbo  Radio  Telescope",  M.D.  Papagiannis  and 
F.L.  Wefer;  Bull.  A AS ,  9,  617,  1977. 

"Observations  of  Enhanced  Radio  Emission  at  15  GHz  from 
a  Coronal  Hole  Region",  F.L.  Wefer  and  M.D. 
Papagiannis;  Bull.  AAS,  9,  617,  1977. 

"The  Radio  Spectrum  of  Coronal  Holes",  H.D. 
Papagiannnis  and  F.L.  Wefer;  EOS,  5_£ ,  3  6  7  ,  1  973  . 


Synoptic  Charts  of  Coronal  Soles  and  9.1  cm  Data  from 
'lav  Through  August  1971",  K . L .  Uefer,  .  1.;. 

i’a  !U'>  i  a  on  i  s  ,  :l.L.  Van  Steenbcn>,  and  T.'l.  Varner; 

Hull.  A AS ,  H),  6  8  A ,  1  978  . 

Solar  Wind  Streams  at  High  Solar  Latitudes  and  Their 
Correlation  to  Coronal  Holes",  M.D.  Papagiannis 
and  K.8.  Baker,  EOS,  h_l_,  1  095  ,  1980  . 


-It  - 
Chapter  1 

RADIO  OBSERVATIONS  OF  CORONAL  HOLES  FROM  THE  SKYLAB  PERIOD 

1.1.  The  Radio  Spectrum  o f  Coronal  Hole 

Although  by  1977  a  number  of  groups  had  observed  coronal 
holes  at  radio  wavelengths  (Du  Ik  and  Sheridan,  1974;  Furst 
and  llirth,  19  7  5  ;  Wefer,  et  al.  ,  1  976  ;  Wefer  and  Blelweiss, 
1976;  Kundu  and  Liu,  1976;  Chiuderi-Drago,  1974;  Lantos 
and  Avingnon,  1975;  Covington,  1976),  no  one  had  analysed  a 
single  coronal  hole  at  many  different  wavelengths  covering  a 
wide  range  of  the  radio  spectrum.  During  the  eight  day 
period  from  June  25,  1973  to  July  2,  1973  the  coronal  hole 
designated  by  the  Skylab  Workshop  on  coronal  holes  as  CHI 
was  visible  on  the  solar  disk.  On  June  28,  1973  solar  radio 
maps  were  made  independently  by  several  groups  at  several 
different  radio  observatories  and  different  radio 
wavelengths:  3.6  inm  and  2.0  cm  -  La  Posta  Observatory,  3.8 
cm  -  NS ROC  Haystack  Observatory,  9.1  cm  -  Stanford  Radio 
Astronomy  Institute,  21  cm  -  University  of  Sydney  La  Fleurs 
Radio  Astronomy  Field  Station.  On  subsequent  passages  of  the 
coronal  hole,  it  was  observed  at  three  additional 
wavelengths:  2.8  cm  -  Bonn  100-m  radio  telescope,  1.8  m  and 
3.75  m  -  Culgoora  r ad i ohe 1 i og r ap h . 

Wefer  and  Papagiannis  (1977)  analysed  these  radio  maps 
in  detail  and  presented  a  radio  spectrum  of  this  hole  by 
choosing  a  point  inside  the  coronal  hole  and  normalizing  the 
brightness  temperatures  at  each  radio  wavelength  with  the 
corresponding  brightness  temperatures  of  the  solar  disk 
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outside  the  coronal  hole  at  the  same  radial  distance  from 
the  center  of  the  disk.  The  combined  observations  were  used 
to  produce  the  radio  spectrum  shown  in  Figure  1.  This 
analvsis  indicated  an  enhancement  of  approximately  2%  in  the 
radio  emission  around  2.0  cm,  which  was  consistent  with  the 
previous  observations  by  Wefer,  et  a  1  .  (1976).  At  2.8  cm  the 
brightness  temperature  of  the  coronal  hole  is  beginning  to 
be  clearly  less  than  the  solar  disk  value,  and  by  21  cm  the 
brightness  temperature  has  decreased  to  about  30%  below  the 
normal  disk  value.  The  minimum  of  the  spectrum,  i.e.  the 
lowest  value  of  the  brightness  temperature  inside  the 
coronal  hole  relative  to  the  brightness  temperature  outside 
the  coronal  hole,  appears  to  occur  somewhere  between  21  cm 
and  1  m . 

The  small  enhancement  (about  2%)  at  2.0  cm  was  then  and 
still  remains  a  controversial  point.  The  enhancement  appears 
to  be  quite  definite  in  the  1973  observations.  The  question 
of  a  possible  enhancement  at  2.0  cm  in  other  coronal  holes 
will  be  discussed  in  more  detail  in  Chapter  A  when  we 
consider  the  1978  observations. 

One  other  important  feature  of  our  analysis  of  CHI  radio 
data  was  the  discovery  that  the  coronal  hole  could  indeed  be 
observed  in  the  9.1  cm  radio  maps  from  Stanford.  The 
Stanford  people  themselves  had  claimed  to  be  unable  to  see 
coronal  holes  In  their  data.  The  fact  that  our  radio 
spectrum  indicates  that  one  ought  to  be  able  to  see  coronal 
holes  around  9  cm  as  10-20%  depressions  in  the  brightness 
temperature  relative  to  the  solar  disk  background  lead  us  to 
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go  brick,  through  all  the  9.1  c  m  Stanford  naps  that  were 
available  during  the  Skylab  period  in  a  search  for  coronal 
holes.  The  results  of  that  search  are  presented  in  the  next 
sect  ion . 

A  detailed  description  of  our  work  on  CHI  is  available 
in  our  Scientific  Report  No.  1,  A FG L- TR- 7 7-0 29  2  ,  Dec.  1977 
by  K.L.  Uefer  and  T.D.  Papagiannis. 


1.2.  Identification  o  f  Coronal  doles  on  the  Stanford  9 . 1  cm 
Synoptic  Charts 

The  period  from  !tay  28,  1  9  7  3  to  August  1  1  ,  1  973  was  an 
ideal  period  for  the  study  of  coronal  holes  because  it  was 
the  only  time  at  which  solar  X-Ray  data  from  ATH-Skylab  and 
9.1  cm  radio  maps  from  the  Stanford  Radio  Astronomy 
Institute  were  simultaneously  available.  In  addition,  it  was 
a  period  near  solar  minimum  when  large  transequatorial 
coronal  (voles  were  present  on  the  sun. 

Our  study  of  the  radio  spectrum  of  coronal  hole  Cdl 
(b'efer  and  Papagiannis,  1977)  had  demonstrated  that  a 
coronal  hole  could  easily  be  seen  at  9.1  cm.  Graf  and 
liracewell  (1975)  had  produced  synoptic  charts  of  the  9.1  cm 
data  for  each  solar  rotation  from  June,  1962  up  to  August, 
1973,  when  the  9.1  cm  system  was  shut  down.  In  their  charts 
the  emphasis  was  placed  on  active  regions  and  the  lowest 
contoured  level  was  at  a  brightness  temperature  of  40,000  K. 
The  central  disk  brightness  temperature  for  the  quiet  sun  at 


9.1  cm  Is  approximately  30,000  K ,  and  since  at  this 
wavelength  coronal  holes  manifest  themselves  as  10-20X 
depressions  in  the  brightness  temperature  relative  to  the 
quiet  sun  background,  it  is  understandable  why  the  synoptic 
charts  produced  by  Graf  and  'Iracewell  did  not  show  any 
coronal  holes  . 

Graf  and  Bracewell  kindly  provided  us  with  the  computer 
tapes  of  the  daily  data  of  their  9.1  cm  observations,  and 
from  these  data,  we  produced  our  own  synoptic  charts, 
emphasizing  contours  near  and  below  the  quiet  sun  level.  The 
basic  procedure  was  a  simple  one.  On  each  day,  we  extracted 
the  data  corresponding  to  the  central  meridian  and  collected 
30  days  worth  of  data  to  form  a  grid  from  which  a  synoptic 
contour  map  was  drawn  covering  slightly  more  than  one 
Carrington  rotation.  Since  this  method  produced  charts  based 
only  on  central  meridian  data,  they  are  suitable  only  for 
studying  long-lived  coronal  holes.  As  a  check  on  the  results 
obtained  in  this  fashion,  we  also  constructed  a  few  charts 
using  data  just  to  the  left  or  right  of  the  central 
meridian.  These  charts  showed  the  same  general  features  as 
the  central  meridian  charts,  but  differed  somewhat  in  the 
finer  (and  shorter-lived)  details. 

A  typical  synoptic  chart  is  shown  in  Figure  2.  The 
shaded  areas  show  the  positions  of  coronal  holes  as 
determined  from  the  ATM-Skylab  data  (Nolte,  et  al.,  1976b) 
for  rotation  1602.  In  this  chart,  CHI  (near  longitude  zero) 
and  CHS  (near  longitude  260°)  are  clearly  manifested  as 
depressions  in  the  9.1  cm  emission.  Coronal  hole  CH3 ,  as 
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well  as  several  smaller  coronal  holes,  did  not  show  up  in 
the  9.1  cm  chart,  probably  because  their  sizes  .ire  smaller 
than  the  beamwidth  of  the  Stanford  antenna  and  the 
additional  side-lobe  interference  from  nearby  active 
r eg i ons  . 

A  complete  description  of  this  work  is  included  in  our 
Scientific  Report  No.  2,  AFGL-TR- 78-0284  ,  January  1979  by 
F.L.  Wefer,  M.D.  Papagiannis,  T.  Varner  and  M .  Van 
Steenberg,  and  in  the  paper,  “Synoptic  Charts  of  Solar  9.1 
cm  and  Coronal  Hole  Data"  by  F.L.  Wefer  and  M.D.  Papagiannis 
(1980). 
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Chapter  2 

H1CH  LATITUDE  CORONAL  HOLES  AND  T'lEIR  CORRELATION  TO 
HIGH  LATITUDE  INTERPLANETARY  SCINTILLATION 

2.1.  Tli  e  Correlation  o  f  IPS  and  Coronal  Hole  Data 

Although  the  correlation  between  large  equatorial 
coronal  holes  and  high  speed  solar  wind  streams  observed 
near  the  earth  has  been  well  established  (Krieger,  et  al., 
1971;  Kopo  and  Orrall,  1  9  7  7  ;  llundhausen,  1  97  7  ),  the 
generalization  that  all  coronal  holes  are  sources  of  high 
speed  solar  wind  streams  is  not  well  established.  In  fact, 

•  lundhausen  (  1  9  7  7  )  has  suggested  that  the  correlation  may  not 
hold  true  for  smaller  coronal  holes. 

A  general  correlation  between  large  coronal  holes  and 
high  solar  wind  velocities,  derived  from  interplanetary 
scintillation  (IPS)  measurements,  was  obtained  in  the  70's 
at  the  University  of  California,  San  Diego  by  Kickett,  et 
al.  (1976).  They  found  that  about  half  of  the  large  coronal 
hole  passages  seen  over  several  solar  rotations  were 
associated  with  high  speed  solar  wind  streams. 

Recently,  Coles,  et  al.  (1980)  have  correlated  the 
average  solar  wind  speed  observed  at  high  solar  latitudes 
with  the  size  of  the  polar  coronal  hole.  They  have  shown 
that  as  the  solar  cycle  approaches  maximum,  the  polar 
coronal  hole  decreases  in  size  and  the  average  solar  wind 
speed  at  high  solar  latitudes  diminishes. 

Nevertheless,  a  definitive  correlation  showing  a 
connection  between  specific  high  latitude  coronal  holes  and 


high  speed  solar  wind  streams  corresponding  to  similar  high 
solar  latitudes  had  not  been  attempted.  To  this  end  we  made 
a  more  detailed  correlation  analysis  between  the  velocities 
of  solar  wind  streams  high  above  or  below  the  ecliptic 
obtained  from  IPS  data,  and  the  positions  and  shapes  of  high 
solar  latitude  coronal  holes  obtained  from  coronal  hole 
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synoptic  maps  constructed  from  the  10330  A  Helium 
spect  roheliograias  . 

For  the  past  several  years  there  have  been  only  two 
groups  making  daily  interplanetary  scintillation 
measurements,  namely,  the  group  headed  by  Coles  at  UCSD  and 
the  group  headed  by  Watanabe  and  Kaklnuma  at  the  Toyokawa 
Observatory  In  Japan.  We  utilized  both  sets  of  data  In  our 
correlation  in  order  to  provide  as  complete  a  coverage  in 
time  as  possible  and  whenever  possible  to  provide  a  check  of 
one  set  against  the  other.  Our  study  covered  the  two  and  a 
half  year  period  from  January  1,  1977  to  June  30,1979.  For 
the  first  year  and  a  half  (Jan. I,  1977  to  June  30,1973)  we 
had  IPS  data  from  both  UCSD  and  Toyokawa.  For  the  last  year, 
however,  we  only  had  data  from  UCSD  because  the  Toyokawa 
data  had  not  yet  been  processed  and  published. 

The  position  and  shapes  of  coronal  holes  were  determined 
from  the  published  monthly  synoptic  charts  (Solar 
Geophysical  Data  Prompt  Reports)  made  by  Dr.  J.  Harvey  from 
from  the  Kitt  Peak  Helium  10830  %  spect rohel iograms .  It 
should  be  noted,  however,  that  because  of  the  difficulties 
involved  in  interpreting  these  s pcc t rohe 1 i ogr ams ,  there  is 
some  uncertainty  in  the  details  and  exact  boundaries  of  the 


coronal  holes  shown  in  the  monthly  synoptic  charts.  In 
addition,  for  the  coronal  holes  present  In  July,  1978,  we 
have  also  made  use  of  our  radio  observations  to  determine 
the  structure  and  boundaries  of  the  coronal  holes.  This  work 
is  discussed  in  detail  in  Chapter  4  and  the  radio  maps  are 
shown  in  Figures  23  -  44. 

The  method  of  analysis  was  a  simple  one.  If  the  point  of 
maximum  scintillation,  P  ,  for  a  particular  radio  source  (see 
Figure  3)  corresponded  to  a  solar  latitude  that  would 
intersect  the  boundaries  of  a  coronal  hole,  we  looked  at  the 
IPS  data  to  determine  if  a  high  speed  solar  wind  stream  was 
present.  We  would  expect  the  coronal  hole  to  produce  a 
scintillation  effect  (a  high  speed  wind  stream)  when  it  is 
pointed  toward  the  point  of  maximum  scintillation.  The 
actual  effect,  however,  would  be  delayed  by  approximately 
one  day,  which  is  the  time  it  would  take  for  the  high  soeei 
wind  stream  to  advance  from  the  coronal  hole  to  the 
scintillation  point.  The  actual  magnitude  of  the  delay  is  a 
function  of  the  solar  wind  speed  and  the  distance  of  the 
scintillation  point  from  the  sun  and  can  be  computed  through 
an  iterative  procedure. 

To  determine  the  correlation  we  had  to  define  a  set  of 
criteria  for  determining  how  good  a  particular  correlation 
was.  We  defined  four  categories  for  the  correlation: 

1)  Very  Good  -  High  speed  solar  wind  streams  were 
observed  on  the  expected  dates.  The  increase  in 
wind  speed  was  significant  (25X  or  more).  If  the 
boundary  of  the  coronal  hole  was  such  that  the 
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scintillation  point  passed  in  and  out  of  the 
coronal  hole  area  more  than  once,  this  should  be 
reflected  in  the  IPS  data. 

2)  Good  -  A  significant  increase  in  the  solar  wind, 
followed  by  a  corresponding  decrease,  must  have 
occurred  within  2  days  of  the  expected  dates. 

3)  Fa 1 r  ~  The  solar  wind  speed  must  have  been  higher 
than  average  during  most  of  the  expected  period 
(with  a  tolerance  of  2  days). 

4)  Poor  -  The  solar  wind  speed  stayed  near  the 
average  value  or  decreased  during  the  dates  of 
expected  high  spjeed  solar  wind. 

Utilizing  these  criteria,  we  were  able  to  analyse  a  total  of 
26  cases  and  classify  the®  in  one  of  the  above  four 
categories.  We  were  unable  to  classify  seven  other  cases  due 
to  a  lack  of  all  the  needed  data.  As  can  be  seen  from  the 
histogram  of  the  26  cases  analysed  (Figure  4),  the 
correlation  in  about  502  of  the  cases  was  either  "Very  Good" 
or  "Good",  while  in  fewer  than  25%  of  the  cases  the 
correlation  was  "Poor". 

Considering  the  assumptions  (radial  propagation  of 
solar  wind  streams,  constant  solar  wind  velocity  along  the 
radial  path,  etc.)  involved  in  this  analysis,  the  sizable 
uncertainties  in  the  solar  wind  speeds  deduced  from  IPS 
measurements,  as  well  as  the  uncertainties  in  the  positions 
and  boundaries  of  coronal  holes  deduced  from  the  10830  A 
synoptic  charts,  the  overall  correlation  is  surprisingly 
good.  We  conclude  that  it  is  certainly  true  that  a 
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significant  fraction  (probably  most)  of  the  high 
latitude  coronal  holes  are  sources  of  high  speed  solar  wind 
streams.  Nevertheless,  it  is  also  quite  possible  that,  like 
the  equatorial  coronal  holes,  some  high  latitude  coronal 
holes  do  not  produce  high  speed  solar  wind  streams  on  a 
continuous  basis,  although  at  various  stages  in  their 
evolution,  they  may  do  so.  It  is  also  quite  evident  that  the 
strength  of  the  solar  wind  streams  is  related  to  the  size  of 
the  coronal  hole.  It  is  quite  possible  therefore  tit  at  some 
of  the  smaller  streams  are  totally  missed,  either  because  of 
their  relatively  short  duration  and  our  limited  means  of 
observation,  or  because  they  get  dissipated  or  diverted  in 
interplanetary  space  before  they  reach  the  expected  point  P 
of  maximum  scintillation. 

A  complete  description  of  this  study,  along  with  the 
data  used  is  given  in  our  Scientific  Report  No.  3, 
APGL-TR-30-0123,  April,  1980  by  K .  15 .  Baker  and  M.D. 
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Chapter  3 


OBSERVATIONS  OF  CORONAL  HOLES  AT  ARECIBO  AND  HAYSTACK 

IN  SEPTEMBER  OF  1977 

3.1.  Coronal  Hoi e  Observations  a  t  11.5  and  2 1  cm  with  the 
A r e c 1 b o  Radio  Telescope  1 n  19  7  7 

During  the  five  day  period  from  August  31  to  September 
4,  1977,  we  obtained  seven  maps  of  the  sun  at  11.5  cm  and 
21  cm  using  the  Arecibo  radio  telescope.  These  solar  maps 
were  made  by  two  different  methods.  In  the  drift  scan  method 
the  telescope  remained  stationary  and  the  sun  was  allowed  to 
drift  past  the  beam  of  the  antenna.  Each  scan  took  about  4 
minutes  and  the  total  map,  consisting  of  21  scans,  took 
about  1.5  hours.  In  the  second  method  the  antenna  was  driven 
at  twice  the  siderial  rate  while  executing  a  boustrop hedonic 
raster.  This  method  made  it  possible  to  cut  the  time 
required  for  each  map  in  half  and  allowed  us  to  make  naps  at 
both  11.5  cm  and  21  cm  during  the  approximately  two  hour 
period  around  local  noon  that  the  sun  can  be  observed  with 
the  Arecibo  radio  telescope  at  this  time  of  the  year.  We 
were  thus  able  to  make  seven  solar  maps  during  this  five  day 
observing  period. 

This  was  the  first  time  that  the  Arecibo  radio 
telescope  was  used  to  make  solar  maps.  The  Arecibo  dish  was 
of  unique  value  for  this  project  because  of  Its  high 
sensitivity  (“1%  of  solar  disk  background)  and  its  high 
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angular  resolution  (3—4  arcnln)  at  decimeter  (11.5  and  21 
cij)  wavelengths.  The  discovery  that  coronal  holes  becone 
clearly  visible  at  these  wavelengths  confirmed  the 
importance  of  the  Arecibo  radio  telescope  for  decimetric 
solar  observations. 

The  results  of  this  project  are  shown  in  Figure  5,  which 
shows  five  of  the  solar  maps  made  during  this  five  day 
perioJ.  The  coronal  hole  boundaries  indicated  on  each  map 

o 

were  obtained  from  the  Kitt  Peak  10830  A  spectroheliograms, 
which  were  provided  by  Dr.  J. Harvey  of  the  Kitt  Peak 
Observatory.  These  spectroheliograms  are  used  for  an  optical 
identification  of  coronal  holes  in  the  absence  of  X-ray  or 
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UV  observations.  The  10330  A  spect rohel i ogram  for  September 
2,  1^77  is  also  shown  in  Figure  5. 

The  coronal  hole  in  the  southern  hemisphere  is  very 
obvious  and  is  very  similar  in  shape  and  position  to  the 
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coronal  hole  seen  in  the  10830  A  s pe c t rob e 1 i og r am .  The 
coronal  hole  In  the  northern  hemisphere  was  much  narrower 
anti  was  therefore  more  difficult  to  observe  at  radio 
wavelengths.  Nevertheless,  a  depression  in  radio  emission  is 
clearly  visible  in  the  region  of  the  north  pole. 
Unfortunately,  the  picture  in  the  north  is  confused  by  the 
presence  of  several  filaments  (indicated  by  dark  wiggly 
lines)  which  also  show  up  as  depressions  at  radio 
wavelengths  (Straka,  Papagiannis  and  Kogut,  1975). 

This  work  is  more  fully  reported  in  the  article, 
"Coronal  holes  at  11.5  and  21  cm  observed  with  the  Arecibo 
radio  telescope",  In  Nature  by  M.D.  Papagiannis  and  F  L. 
W  e  f  e  r  ,  l  9  7  .3  . 


1.2.  C  o  r  o  n  a  1  'do_l_e  Obse  rv.it  ions  .1 1  II  a  yst.ick  a  n  d  L  a  P  o  s  t  a  1  n 
197  7 

In  addition  to  the  Areclbo  radio  naps,  we  were  also  able 
during  the  same  period  to  make  radio  maps  at  3.81  cm  and 
1.94  cm  using  the  Haystack  Observatory  radio  telescope. 
These  maps,  a  Ion'.;  with  maps  at  2.0  cm  and  0.36  cm  from  the 
La  Posta  Observatory  are  shown  in  Figure  6.  The  position  of 
the  coronal  holes  and  filaments  are  also  indicated  on  these 
maps.  The  coronal  holes  appear  as  slight  depressions  at  3.31 
cm.  At  2.0  cm  and  1.94  cm,  the  holes  anpear  as  slight 
enhancements.  This  apparent  enhancement  will  be  discussed  in 
more  detail  in  Chapters  4  and  6. 

3.3.  The  Rad io  Spec  t  rum  o  f  the  September  1  977  Coronal  Hole 

From  our  observations  at  Areclbo  and  Haystack,  combined 

with  additional  solar  radio  mans  made  by  the  La  Posta 
observatory  in  California,  we  were  able  to  determine  a 
partial  radio  spectrum  of  this  coronal  hole  (Figure  7).  The 
results  are  quite  consistent  with  those  for  CHI  in  1973 
(Figure  1  in  Chapter  1)  particularly  with  respect  to  the 
approximately  2%  enhancement  around  2  cm. 


TENNA  TEI 


FIGURE  7 


Chapter  4 


CORONAL  HOLE  RADIO  OBSERVATIONS  DURINC  JULY  1973 

4.1.  International  Effort 

From  the  previous  work  it  had  become  clear  that 
simultaneous  observations  of  coronal  holes  over  a  wide  r.inj’e 
of  radio  wavelengths  would  allow  us  to  determine  the  radio 
spectrum  of  coronal  holes,  which  would  be  of  great  value  in 
our  effort  to  understand  the  internal  structure  and  the 
physics  of  coronal  holes.  It  would  also  he  useful  to  have  a 
certain  amount  of  redundancy  in  the  observations  to  insure  a 
higher  level  of  dependability  on  the  data. 

Starting  early  in  1978  we  began  to  organize  an 
international  co-operative  effort  to  observe  coronal  holes 
at  many  different  radio  wavelengths.  The  project  was  well 
received  by  many  groups  around  the  world  and  the  tine  period 
selected  for  this  combined  effort  was  July  of  1973.  The 
month  of  July  was  chosen  primarily  because  we  wanted  to 
Include  observations  at  21  cm  and  11.5  cm  from  the  Arecibo 
observatory  which  can  observe  the  sun  only  during  the  summer 
and  allows  the  longest  observing  time  near  the  summer 
solstice.  In  addition,  most  of  the  participating, 
observatories  were  located  in  the  northern  hemisphere  and 
would  be  expected  to  have  better  observing  conditions  during 
the  summer  months.  In  order  to  increase  the  probability  that 
a  coronal  hole  would  be  present  on  the  solar  disk  during  our 
observations,  wc  decided  to  use  two  separate  observing 
periods  of  about  five  days  each  separated  by  about  ten  days. 
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The  two  periods  were  JuLv 

Observatory  Location 

6-10  . 
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these 

ma  p  s 
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summarized  in  Table  1.  Unfortunately  the  radio  maps  from  the 
Toy  ok aw a  observatory  and  the  University  of  British  Columbia 
had  severe  sidelobe  problems  and  difficulties  with  the 
processing  of  their  radio  data  which  made  these  maps 
virtually  useless  for  studying  coronal  holes.  These  maps, 
therefore,  were  not  Included  In  our  analysis. 

4.2.  First  Period  Observat  ions  (July  6-  1  0,  1978) 

The  first  period  chosen  for  the  radio  observations  was 
from  July  6  to  July  10,  1978.  Not  all  observatories  were 
able  to  make  observations  for  the  entire  period,  but  there 
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is  nonetheless  a  reasonable  overlap  of  dates  from  all 
observatories.  Table  2  1 tsts  the  solar  parameters,  P,  B0  , 
and  L0  (Carrington  Longitude  of  the  Central  Meridian)  at 
0000  UT  for  each  of  the  six  days  of  the  first  period. 


TABLE  2 

Solar  Parameters  at  0000  UT 


DATE 

ha 

io 

P 

7-6 

23  5 

+  3.4 

-0.6 

7-7 

272 

3 . 5 

-0.1 

7-3 

259 

3.6 

+  0.3 

7-9 

246 

3.7 

+0.3 

7-10 

232 

3.8 

+  1  .2 

Table  3  lists  the  participating  observatories  and  indicates 
the  wavelengths  and  the  days  on  which  each  observatory  made 
its  observa  t ions . 


TABLE  3 

Participating  Observatories  -  First  Period  Observations 


July  Date 


Observa  tory 

Wavelength(cm)  6 

7 

8 

9 

10 

La  Pos  t  a 

0.86 

X 

X 

X 

X 

X 

La  P o s  t  a 

2.00 

X 

X 

X 

X 

'lav  stack 

0.69 

X 

X 

day  stack 

1  .38 

X 

X 

X 

Kay  stack 

1.93 

X 

X 

X 

Haystack 

3.81 

X 

X 

1 1  a  pe  t i nga 

0.7  0 

X 

X 

X 

I tapet inga 

1  .39 

X 

X 

X 

A  r  e  c  i  b  o 

11.5 

X 

X 

X 

A  r e  c i bo 

2  1  .4 

X 

X 

C  u  1  g  o  o  r  a 

137 

X 

Culgoora 

375 

X 

X 

As  is  clear 

f  rom  Table  3 

,  July  8 

is  the 

date 

on 

which 

mos  t  comple  te 

coverage  is 

available 

• 

The  first  period  turned  out  to  be  a  poor  one  for  the 
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purpose  of  observing  coronal  holes.  The  Helium  10830  A  data 
show  only  the  usual  small  polar  depressions.  The  helium 
synoptic  map  from  Kitt  Peak  for  half  a  Carrington  rotation 
centered  approximately  on  July  7  is  shown  in  Figure  8.  There 
are  no  coronal  holes  but  there  are  several  active  regions 
present.  These  active  regions  can  also  be  clearly  seen  in 
the  Ha  and  Calcium  pictures  of  the  sun  for  this  period.  The 
ila  picture  from  Sacramento  Peak  and  the  Ma  cMa  t  h- II  u  1  be  r  t 
Calcium  Report  for  July  8  are  also  shown  in  Figure  3.  The 
three  largest  active  regions  are  McMath  15335,  15391,  and 
15393.  The  Ha  picture  shows  a  very  clearly  defined  filament 
in  the  northern  hemisphere  at  an  approximate  Carrington 
longitude  of  240°  .  Thus,  during  the  first  period  this 
filament  is  always  on  the  east  side  of  the  central  meridian. 
It  is  important  to  note  the  position  of  this  filament, 
since,  as  mentioned  above,  large  filaments  appear  very 
similar  to  coronal  holes  on  radio  maps. 


Radio  Maps  for  July  6-10, 1978 

Figures  9  -  20  show  representative  maps  from  the  first 
period  observations.  Although  some  details  vary  from  day  to 
day,  the  important  features  are  relatively  constant  over  the 
five  day  period  from  July  6  to  July  10. 

The  shortest  wavelength  data  available  are  the  0.69  cm 
maps  from  the  Haystack  observatory.  At  this  very  short 
wavelength  all  the  information  comes  from  deep  in  the 
chromosphere  and  all  solar  features,  including  the  active 
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regions,  show  very  Little  contrast  with  the  stir  round  i  ns; 
quiet  area.  Still,  the  active  regions  can  be  seen,  as  well 
as  the  suggestion  of  a  slight  depression  corresponding  to 
the  position  of  the  filament  in  the  northeast.  The 
Itapetinga  maps  at  0.70  cm  for  the  first  period  were  marred 
bv  spurious  readings.  We  have  cleaned  up  the  map  by 
interpolation  to  remove  the  spurious  data,  but  the  resultant 
maps  show  few  features. 

The  'laystack  l  .Vi  cm  maps  and  the  Itapetinga  1.39  cm 
naps  noth  clearly  show  the  active  regions  and  the  filament 
in  the  northeast.  The  Itapetinga  maps  were  again  marred  by 
spurious  data,  but  in  regions  where  we  can  be  sure  the 
Itapetinga  naps  are  reasonably  clean,  tbe  features  are  quite 
consistent  with  the  Maystack  maps. 

The  La  Post  a  2.0  cm  maps  and  the  Haystack  1.93  cm  maps 
are  also  quite  consistent.  Hoth  show  the  active  regions  as 
strong  enhancements  and  the  filament  as  a  slight  depression. 
Hie  la /stack  maps  also  have  a  suggestion  of  a  depression 
slightly  east  of  the  central  meridian  in  the  far  south.  The 
Sac  Peas  ’ip  map  indicates  the  presence  of  a  small  filament 
l n  this  region,  which  Is  probably  the  source  of  this 
depress  Ion . 

The  Arecibo  maps  at  11.5  cm  and  21.4  cm  show  the  active 
regions  and  the  filament  in  the  northeast  quite  clearly. 

The  long  wavelength  maps  from  the  Culgoora  Observatory 
(  1  . -i  7  m  and  3.75  m)  show  much  less  detail  than  the  other 
radio  -naps.  The  active  regions  can  be  seen  at  both 
wavelengths,  but  at  the  shorter  wavelength  the  maps  are 
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confused  by  negative  sidelobes  from  these  active  regions.  \n 
additional  complication  in  interpreting  these  maps  is  that 
the  sun  was  very  low  in  the  skv  during  these  July 
observations  at  Culgoora.  Refraction  of  these  long 
wavelength  radio  waves  through  the  earth's  atmosphere 
displaces  the  radio  image  of  the  sun  from  the  epheraerls 
position,  making  it  difficult  to  determine  the  precise 
position  of  the  observed  features  on  the  solar  disk. 
Unfortunately  these  maps  were  only  available  in  the  form  of 
finished  contour  maps  and  not  In  the  form  of  data  points. 
This  made  it  impossible  for  us  to  re-contour  the  maps  with 
emphasis  on  contours  near  the  quiet  sun  level. 

4.3.  Second  Period  Observations  (July  15-13,  1978) 

The  coverage  of  the  radio  spectrum  during  the  second 
period  was  more  complete  than  for  the  first  period.  Table  4 
lists  the  solar  parameters  for  the  four  days  of  the  second 
period  and  Table  5  lists  the  participating  observatories, 
wavelengths,  and  the  days  each  observatory  made  Its 
observations.  July  16  and  17  are  the  dates  with  the  best 
cove  rage . 

TABLE  4 

Solar  Parameters  at  0000  UT 


DATE 

k  0 

lo 

P 

7-15 

166 

4.3 

3.4 

7-16 

1  53 

4.4 

3.9 

7-17 

140 

4.5 

4.4 

7-18 

127 

4.6 

4.3 

The  Helium  10330  A  s pec t r oh e 1 l ogr am  for  July  16  is  shown 
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in  Figure  21,  and  the  Helium  synoptic  clinrt  for  the  half 
C.irri,i",t  ion  rotation  centered  approximately  on  July  16  is 
shown  in  Figure  22.  The  outlines  of  the  coronal  holes  are 
also  shown  as  they  were  drawn  by  Jack  Harvey,  of  the  Kitt 
Peak  Observatory,  who  kindly  provided  us  with  these  maps. 
The  helium  data  show  a  large  coronal  hole  extending  from  the 
south  pole  toward  the  equator,  a  small  hole  at  latitude  1 0° S 
and  a  somewhat  larger  hole  at  latitude  40°N.  \  number  of 
active  regions  are  also  present. 


T  A  3  L  f.  5 


Participatin',  Observatories  -  Second  Period  Observations 


Observatory 

Wave  length(cm) 

l  5 

July 

16 

Date 

17 

13 

La  Post  a 

0.8  6 

X 

X 

X 

X 

L  a  Post  a 

2.00 

X 

X 

X 

X 

1  a  y  s  t  a  c  k 

0.69 

X 

X 

X 

Haystack 

1.38 

X 

X 

X 

i  a  y  s  t  a  e  k 

1.93 

X 

X 

X 

1 1  a  y  s  t  a  c  k 

3.8  1 

X 

X 

X 

I  t  a  p  e  t in  g a 

0.70 

X 

X 

X 

l  t.inet  i  n g a 

l  .39 

X 

X 

X 

A  r  e  c  i  b  o 

11.5 

X 

X 

X 

A r e  c i ho 

Ju  1  !•  o  o  r  a 

2  1.4 

137 

X 

X 

X 

V 

X 

C  u  1  >o  ora 

37  5 

X 

X 

X 

The  Sac  Peak  da 

picture 

(Figure 

22) 

for 

July 

16  also 

hows  several  strong 

act  l  v  e 

regions 

a  s 

well 

as 

several 

filaments.  Also  shown  in  Figure  22  is  the  McMa t h-Uul be rt 
Calcium  deport.  The  most  prominent  active  region  is  flcMath 
1  540  1  . 

The  coronal  hole  in  the  south  lies  in  a  relatively 
quiet  area  of  the  sun  which  should  make  it  easily  detectable 


-  51  - 

at  radio  wavelengths.  On  the  other  hand,  the  coronal  hole  In 
the  northern  hemisphere  lies  very  close  to  the  large 
filament  between  Me ‘lath  15410  and  McMath  15415,  which  makes 

It  difficult  to  identify  the  coronal  hole  at  radio 
wave le  ng t  hs  . 


i 


i 


i 
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H at!  i  o  :1a p s  for  .Inly  1  r>  -  1  8  ,  1978 

Figures  23  -  44  show  some  of  the  naps  for  the  second 

period  observations.  Special  emphasis  has  been  placed  on  the 
maps  for  the  two  days  July  16  and  July  17,  since  these  two 
Jays  represent  the  days  for  which  we  have  the  most  complete 

wavelength  coverage. 

As  with  the  first  period  maps,  the  shortest  wavelength 
maps  show  little  contrast.  The  Haystack  0.69  cm  maps  are 
virtually  devoid  of  features  other  than  the  strong  active 
regions.  The  Itapetinga  maps  at  0.70  cm  show  a  distinct 
depression  in  the  southern  hemisphere  corresponding  in 
position  to  the  coronal  hole.  In  contrast,  however,  the  0.36 
cm  maps  from  La  l’osta  show  only  the  filament  and  coronal 
hole  in  the  northeast  as  depressions. 

"lie  Haystack  maps  at  1.33  cm  also  show  distinct 
differences  from  the  Itapetinga  maps  at  1.39  cm.  The 
Haystack  maps  show  clear  depressions  in  both  the  north  and 
the  south.  In  both  cases,  however,  the  position  of  the 
depression  is  shifted  to  the  west  of  the  position  of  the 
coronal  hole  as  determined  from  the  Helium  10330  X  charts. 
This  apparent  shift  In  the  position  had  also  been  noted  in 
previous  observations  (We  for  and  Papagiannis,  1977). 

The  Haystack  maps  also  show  a  depression  on  the  east  limb 


in  the 

southern  hemisphere. 

The 

Ha 

pictures  show 

seve  ra 1 

small 

filaments  in  this  regi 

on  of 

the  east  limb  which 

a  re 

probably  unresolved  in  the 

r  a  d  i 

O 

maps  and  show 

up  as 

a 

single 

extended  depression. 

The 

Itapetinga  maps  at 

1  .39 

c  m 

show 

only  a  slight  suggestion 

o  f 

the  depression 

i  n 

t  he 
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north,  while  the  depression  in  the  south  appears  to  he 
merely  an  extension  of  the  large  depression  coining  from  the 
east  limb.  Presumably  the  differences  between  the  Haystack, 
and  Itapetinga  maps  reflects  the  difference  in  angular 
resolution  of  the  two  antennas.  At  Haystack  the  resolution 
is  about  1.6  arcrain  at  1.38  cm  while  at  Itapetinga  the 
resolution  is  only  about  4  arcmln.  Thus,  the  nearby  active 
regions  mask  the  coronal  holes  in  the  Itapetinga  maps, 
whereas  the  Haystack  antenna  was  able  to  resolve  these  as 
separate  features. 

The  maps  at  1.93  c.i  (Haystack)  and  2.0  c.i  (La  Posta) 
uoth  show  slight  depressions  corresponding  to  the  coronal 
hole  in  the  south.  As  with  the  1.38  cm  maps,  the  position  of 
the  coronal  hole  appears  to  be  shifted  to  the  west  of  its 
position  in  the  Helium  10830  X  maps.  Indeed,  this  westward 
shift  is  a  consistent  feature  of  the  longer  wavelength  maps. 
In  the  north,  the  Haystack  maps  show  the  coronal  hole  as  a 
distinct  depression.  Slightly  west  of  the  depression, 
however,  is  a  small  region  of  enhanced  emission.  This 
enhancement  does  not  correspond  to  any  clear  features  In  the 
Ha  or  Calcium  s pe c t r ohe 1 i og r ams ,  but  it  does  appear  to  be 
real . 

The  La  Posta  2.0  cm  maps  indicate  a  slight  enhancement 
In  the  region  of  the  northern  coronal  hole.  This 
disagreement  between  Haystack  maps  and  La  Posta  maps  may  be 
the  result  of  differences  In  angular  resolution.  The  regions 
of  enhanced  emission  seen  In  the  Haystack  maps  as  separate 
features  are  probably  masking  the  slight  depression  of  the 


coronal  hole  in  the  La  1‘osta  maps.  If  slight  enhancements 
near  the  boundary  of  a  coronal  hole  are  common  features, 
this  may  explain  the  previous  results  of  Wefer,  et  al. 
(197b),  Wefer  and  Hleiweiss  (1976),  and  our  own  results  from 
1977  (see  Chapter  3),  where  coronal  holes  appeared  as 
e  nli  a  nc  e  me  n  t  s  of  about  1  -  2  &  at  2.0  cm. 

The  Haystack  maps  at  3.31  cm  show  the  usual  active 
regions  as  well  as  the  coronal  holes  in  both  the  north  and 
the  south.  As  with  the  1.93  cm  maps,  these  maps  show  a  small 
region  of  enhanced  emission  slightly  to  the  west  of  the 
northern  coronal  hole. 

The  Arecibo  maps  at  both  21.4  cm  and  11.5  cm  show  both 
coronal  holes  quite  clearly,  including  the  eastward 
extension  of  the  southern  coronal  hole.  The  depression  on 
the  east  limb  is  also  clearly  visible  as  a  separate  feature. 
As  before,  the  active  regions  are  also  prominent  features. 

As  was  the  case  during  the  first  period  observations, 
the  Culgoora  maps  at  1.37  m  and  3.75  in  are  more  difficult  to 
interpret  than  the  other  maps,  due  to  the  uncertainty  about 
the  position  of  the  solar  disk  and  the  coarse  contour 
intervals  in  the  solar  maps  provided.  At  3.75  m,  the  most 
prominent  feature  is  the  large  brig  lit  region  near  the 
center,  which  is  probably  the  strong  active  region  McMath 
15403.  South  of  this  feature  is  a  slight  depression  which 
corresponds  in  position  to  the  coronal  hole.  As  usual,  the 
1.37  m  maps  show  many  more  features  than  the  3.75  m  maps. 
Again,  sidelobes  create  major  difficulties  In  interpreting 
these  .naps.  There  appear  to  be  two  regions  of  depressed 
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emission  in  the  south,  with  the  eastern  one  corresponding 
more  closely  in  position  to  the  coronal  hole. 

4.4.  The  Spectra  o f  the  Second  Period  Coronal  Holes 

Determining  a  spectrum  for  a  coronal  hole  is  a  rather 
difficult  task  when  there  are  many  active  regions  on  the 
solar  disk.  The  spectrum  is  actually  the  ratio  of  the 
brightness  temperature  inside  the  coronal  hole  at  a 
particular  wavelength  to  the  brightness  temperature  that  we 
would  have  observed  in  the  same  region  if  the  coronal  hole 
had  not  been  present.  The  best  approach  for  determining  the 
values  for  the  quiet  sun  is  to  measure  the  brightness 
temperature  in  a  quiet  area  at  the  same  distance  from  the 
center  of  the  solar  disk  and  at  the  same  latitude  as  the 
coronal  hole.  Unfortunately,  when  there  are  many  active 
regions  present,  such  quiet  areas  become  difficult  to  find, 
especially  since  the  relatively  large  beamwldths  of  the 
antennas  increase  the  area  of  the  solar  disk  affected  by  the 
active  regions. 

This  problem  was  especially  acute  during  the  second 
period,  when  two  coronal  holes,  one  in  the  northern  and  one 
in  the  southern  hemisphere,  were  situated  near  the  central 
meridian  at  approximately  the  same  distance  from  the  center 
of  the  solar  disk.  In  addition,  four  active  regions  were 
nearly  symmetrically  arranged  in  the  four  quadrants  of  the 
solar  disk.  This  made  it  unusually  difficult  to  determine 
the  appropriate  quiet  solar  brightness  temperatures,  which 
introduced  some  additional  uncertainty  in  the  determination 


of  the  radio  spectra. 

In  this  section  we  will  determine  the  spectra  of  the  two 
coronal  holes,  the  one  in  the  south  and  the  one  in  the 
north,  on  two  consecutive  days,  July  If)  and  17.  The 
determination  of  the  brightness  temperature  ratio  will  be 
discussed  for  each  map.  The  Toyokawa  maps  will  not  be  used 
because  severe  sidelobe  problems  make  them  completely 
unreliable . 

T h  e  Spectrum  o  t  the  Northern  Coronal  Hole  July  1 6 ,  1978 

The  coronal  hole  in  the  north  was  an  isolated  hole,  not 
connected  to  the  polar  coronal  hole.  Since  this  hole  does 
not  cover  a  wide  range  of  latitudes,  there  is  no  point  in 
trying  to  get  a  spectrum  at  several  different  points  within 
the  hole.  Indeed,  such  an  attempt  would  probably  simply 
result  in  misinformation  because  near  the  boundary  the 
antenna  beams  would  be  seeing  as  much  area  outside  the  area 
of  tlie  coronal  hole  as  Inside.  We  will  therefore  look  for 
the  point  of  maximum  contrast. 

In  determining  the  spectrum  of  a  coronal  hole  when  we 
had  two  maps  at  very  similar  wavelengths  we  used  the  map 
with  the  higher  resolution.  Thus,  wo  have  chosen  the 
Haystack  map  at  1.93  cm  in  preference  to  the  La  Posta  2.0  cm 
map  anil  the  Haystack  maps  at  1.38  cm  and  0.69  cm  in 
preference  to  the  Itapetinga  maps  at  1.39  cm  and  0.70  cm. 

Table  6  lists  the  observatories;  wavelengths;  the 
quiet  sun  brightness  temperatures,  T  ^  in  relative  units;  the 
coronal  hole  brightness  temperatures,  in  the  same  units; 


and  the  ratios  of  these  two  brightness  temperatures  for  the 


northern 

co'ona 1  hole 

on  duly  lb,  1973 

.  The 

Culgoora  ma  i)  is 

difficult 

to  interpret 

and  the 

brightness  temperature  ratio 

at  375  cm 

must  be  treated  with 

c  a  u t ion 

.  The 

values  obtained 

most  likely  represent 

a  lower 

bound,  and  the 

s  pe  c  t  ra 1  value 

at  375  cm 

could  be  considerably 

closer 

to  1  . 

TABLE 

6 

Spectrum  of  the 

Northern 

Coronal 

Hole  - 

July  16 

Observatory 

X(  cm  ) 

T 

q 

hx 

Th/Tq 

Culgoora 

37  5 

356 

3  06 

0.86 

A  re  c l bo 

2  1  .  A 

80 

55 

0.69 

Arecibo 

1  1  .  A 

85 

7  0 

0.82 

Haystack 

3.81 

1  00 

92 

0.9  2 

Haystack 

1.93 

100 

95 

0.95 

Haystack 

1  .38 

99 

97 

0.9S 

La  Posta 

0.86 

36.5 

36 

0.99 

Haystack 

0.69 

8  5 

85 

1  .00 

Of  particular  Interest  are  the  Haystack  observations. 
The  region  In  the  northern  hemisphere  corresponding  most 
closely  to  the  coronal  hole,  as  determined  from  the  10010  % 
Helium  maps  and  the  Areclho  radio  maps,  is  a  region  of 
distinctly  lower  brightness  than  the  quiet  sun  level. 
However,  at  the  northeastern  boundary  of  the  coronal  hole 
region  the  3.81  cm  map  shows  the  suggestion  of  an 
enhancement.  Furthermore,  at  1.93  cm  and  1.33  cm  these 
regions  appear  as  very  distinct  enhancements  approximately 
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2-3  v.  Sri'.’htcr  than  the  expected  q  u  i  e  t  sun  level.  This 
enhancement  is  consistent  uith  that  previously  reported  by 
'defer,  et  a  1  .  (1976)  and  Wefer  and  P  a  p  a  g  i  a  n  n  i  s  (1977)  for 
coronal  holes  around  2  cm.  Since  the  spatial  resolution  of 
the  La  Post.)  telescope  used  by  'defer  et  al.  is  much  poorer 
titan  the  ‘iavstack  telescope,  it  is  conceivable  that  what  was 
being  observed  were  similar  small  enhancements  near  the 
coronal  hole  boundaries.  Although  no  plage  area  can  he  seen 
in  the  neighborhood  of  this  enhancement,  it  is  quite 
possible  that  near  the  boundaries  of  a  coronal  hole,  closed 
magnet ic  field  lines  lead  to  higher  densities  and  hence  in 
some  cases  to  regions  of  slightly  enhanced  emission. 

f  he  S p ectru m  t  he  Northern  Coronal  hole  -  duly  1 7 ,  19  7  3 

'in  duly  17,  we  have  an  additional  map  from  Culgoora, 
Australia,  but  we  are  missing  the  0.69  cm  map  from  the 
Haystack  observatory.  The  Culgoora  160  MHz  map  is 
complicated  by  sidelobe  problems,  which  make  it  impossible 
to  locate  the  coronal  holes  unambiguously  and  to  determine  a 
relative  brightness  temperature  for  them. 

Table  7  presents  the  results  for  July  17.  The  spectrum 
is  quite  consistent  with  the  results  of  duly  16.  The 
contrast  appears  to  he  slightly  greater  at  11.5,  3.81,  and 
1.38  cm,  but  the  error  bars  for  these  measurements  are  about 
2-3%  so  these  differences  are  not  serious.  We  note  again 
that  the  coronal  hole  appears  as  a  depression  at  1.93  cm  and 
again  there  are  two  small  enhancements  (about  2-32)  near  the 
northeastern  boundary  of  the  coronal  hole. 
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figure  45  shows  the  spectra  for  the  northern  coronal 
hole  for  both  days.  The  solid  line  is-  the  spectrum  for  July 
16,  while  the  dashed  line  is  the  spectrum  for  July  17.  The 
error  bars  were  estimated  directly  from  the  maps. 


85 


The  Spectra  o f  the  Southern  Coronal  Hole  -  July  1 6 

The  southern  coronal  hole  consists  of  two  parts.  The 
primary  hole  is  an  extension  of  the  polar  coronal  hole, 
which  reaches  a  latitude  of  approximately  JO  S.  The  other 
part  is  a  small  coronal  hole  northeast  of  the  main  extension 
at  about  10*  S  latitude. 
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This  small  region  lies  very  near  a  filament,  which  makes 
interpreting  the  maps  rather  difficult.  In  order  to  avoid 
this  difficulty,  we  will  concentrate  on  the  main  coronal 
hole.  Since  this  coronal  hole  covers  a  large  range  of 
latitudes,  we  can  gain  some  additional  information  by 
determining  the  coronal  hole  spectrum  not  only  at  the 
position  of  greatest  contrast,  but  also  at  a  larger  aiming 
parameter  near  the  limb. 

Table  8  presents  the  data  for  two  different  positions 
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in  the  coronal  hole.  The  first  is  at  the  region  of  greatest 
contrast  (about  40  S)  and  the  second  is  near  the 
photospheri c  limb.  In  general,  near  the  limb  the  brightness 
contrast  is  much  less.  However,  it  is  interesting  that  the 
Haystack  map  at  0.69  cm  and  the  Itapetinga  map  at  0.70  cm 
suggest  a  slightly  increased  contrast  between  the  coronal 
hole  and  t  tie  normal  brightness  near  the  limb.  The  Haystack 
mu p  it  1.93  cm  shows  very  little  contrast  near  the  limb, 
except  for  two  small  regions  near  the  coronal  hole.  The  one 
on  the  left  shows  a  depression  of  about  2%,  while  the  one  on 
the  right  shows  an  enhancement  of  about  1-2%.  This  slightly 
enhanced  region  (which  is  almost  directly  on  the  central 
meridian)  may  not  lie  within  the  area  of  the  coronal  hole  and 
nay,  therefore  be  a  region  similar  to  the  slight 
'■'I  i  near  the  northern  coronal  hole.  If,  however, 
this  region  is  within  the  coronal  hole,  the  result  would  be 
consistent  with  the  earlier  results  discussed  in  Chapter  3. 

The  Spectra  o  f  the  Southern  Coronal  Hole  July  1 7 

Table  9  gives  the  coronal  hole  spectra  at  the  points  of 
highest  contrast  and  near  the  limb  for  July  17.  The  Culgoora 
160  MHz  map  is  sufficiently  clear  in  the  southern  hemisphere 
(as  opposed  to  the  northern  hemisphere  -  see  above)  to  allow 
us  to  include  It  In  the  determination  of  the  spectrum.  As 
with  the  northern  coronal  hole,  the  r-sults  are  similar 
enougii  to  the  results  for  July  16  to  be  considered 
consistent.  Figure  46  shows  the  spectra  for  the  southern 
coronal  hole  for  both  days  and  both  positions. 
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Chapter  5 

SPECTRA  AND  CORONAL  MOLE  MODELS 
5.1.  The  Spectra  o f  Coronal  Holes 

Over  the  period  covered  by  this  contract  we  have  been 
able  to  determine  radio  spectra  for  four  different  coronal 
holes.  In  one  case,  the  southern  coronal  hole  in  July,  1975, 
we  were  even  able  to  determine  the  spectrum  at  two  different 
positions  inside  the  hole,  one  near  the  polar  limb  and  the 
other  at  the  point  of  highest  contrast  at  about  4  0°  S 
latitude. 

Although  the  individual  spectra  are  discussed  in  the 
previous  chapters,  it  is  instructive  to  compare  all  of  them 
here  to  see  if  there  are  any  significant  differences  among 
the  different  coronal  holes.  Overall,  the  results  are  quite 
consistent,  and  the  differences  between  coronal  holes  seem 
to  be  within  the  error  bars  of  our  observations  and 
analysis.  The  greatest  area  of  disagreement  occurs  around 
2.0  cm.  The  spectrum  for  CHI  of  1971  and  the  spectrum  for 
the  coronal  hole  of  September,  1977  both  show  a  slight 
enhancement  of  approximately  2%  around  t h  i  s  wavelength.  On 
the  other  hand,  neither  of  the  two  coronal  holes  in  the 
July,  1973  data  show  this  enhancement,  although  in  one  case 
there  was  a  region  near  and  possibly  inside  the  coronal  hole 
that  does  show  an  enhancement.  As  was  pointed  out  in  the 
previous  chaptei,  the  Haystack  data  for  July,  1978  do  show  a 
region  of  enhanced  emission  near  the  northern  coronal  hole 
but  not  inside  the  hole,  while  the  southern  cornal  hole 
shows  a  slight  enhancement  near  the  limb  which  may  or  may 


not  he  Inside  the  coron.il  hole 


It  is  interest  ini»  to  note  that  the  data  for  CHI  indicate 
a  depression  of  about  3%  at  2.8  cm  (Furst  and  Mirth,  1975) 
and  a  depression  of  about  2%  at  3.8  cm.  This  rather  odd  dip 
in  the  spectrum  is  suspicious  and  may  indicate  that  the 
errors  in  determining  the  relative  brightness  temperatures 
were  ,» renter  than  estimated.  If  so,  the  data  are  consistent 
with  no  detectable  difference  in  brightness  temperature 
arounJ  2  cm  between  a  coronal  hole  and  the  quiet  solar 
back  g round . 

Finally,  it  must  be  pointed  out  that  the  spectrum  near 
the  limb  for  the  July,  1978  southern  coronal  hole  shows  much 
less  contrast  when  compared  with  the  spectrum  of  the  same 
coronal  hole  on  the  disk.  This,  however,  is  to  be  expected, 
since  near  the  limb  we  do  not  reach  as  deep  into  the  solar 
atmosphere,  where  the  effects  are  more  pronounced,  as  we  do 
near  the  center  of  the  solar  disk. 


9.2.  I o d e 1 i n g  -  Procedures 

'laving  determined  the  radio  spectrum  of  a  given  coronal 
hoie,  we  now  want  to  use  it  to  obtain  the  changes  that  occur 
in  the  electron  density,  temperature,  and  neutral  particle 
density  profiles  inside  the  coronal  hole.  The  brightness 
temperature  at  any  given  wavelength  Is  a  function  of  the 
temperatures,  electron  densities,  and  neutral  particle 
densities  along  the  ray  path  to  tnc  observer.  At  shorter 
wavelengths  the  emission  originates  deeper  inside  the  solar 
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a  titios  pile  t  e ,  while  at  longer  wavelengths  It  originates  higher 
up.  With  observations  available  from  n  . 7  cm  to  175  cm,  i.e. 
over  more  than  9  octaves,  we  are  able  to  gain  information 
from  as  deep  as  the  lower  chromosphere,  through  the 
important  transition  zone  and  all  the  way  to  the  outer 
corona . 

The  process  of  analyzing  the  data  consists  essentially 
of  using  a  ray  tracing  procedure  to  compute  the  brightness 
temperature  of  the  sun  at  a  given  wavelength  and  radial 
distance  (aiming  parameter)  from  the  center  of  the  solar 
disk.  This  is  done  first  for  the  quiet  sun  and  then  for  a 
sequence  of  modified  solar  atmospheres.  The  ratios  of  the 
computed  brightness  temperatures  for  the  modified  and  the 
quiet  solar  atmosphere  at  a  series  of  wavelengths  are  used 
to  construct  the  radio  spectrum  for  a  particular  modified 
model  of  the  solar  atmosphere.  This  is  then  compared  with 
the  actual  radio  spectrum  of  coronal  holes  deduced  from  our 
observations  in  an  effort  to  determine  the  modified  model 
that  most  closely  matches  the  observed  spectral 
characteristics  of  coronal  holes. 

This  procedure  requires  a  standard  atmosphere  for  the 
quiet  sun  which  can  then  be  modified  to  find  the  changes 
inside  the  coronal  holes.  The  model  we  have  used  for  the 
quiet  sun  atmosphere  is  the  Boston  University  Solar  Model 
(Papagiannis  and  Kogut;  1975,  1976)  which  follows  the 
Vernazza,  et  al.  (1973)  model  for  the  chromosphere  up  to  a 
height  of  2400  km  above  the  photosphere,  and  is  then 
connected  using  the  method  described  by  Wlthbroe  (1970)  to  a 
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Dupree-Goldberg  (1967)  model  for  the  transition  zone  and  the 
corona.  The  electron  density,  the  temperature,  T,  and 

the  collision  frequency,  ,  of  the  model  near  the  critical 
transition  zone  are  shown  in  Figure  47,  while  the  complete 
listing  of  the  different  parameters  of  the  model  Is  given  in 
Appendix  1. 

The  brightness  temperature,  T^ ,  is  given  by  the  integral 
along  the  ray  path 
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where  T  is  the  plasma  temperature  and  T  is  the  optical 
depth.  Since  this  is  actually  calculated  in  discrete  steps, 
the  brightness  temperature  is  approximated  by  the  expression 
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The  equation  for  a  ray  trajectory  is  (Pawsey  and  Smerd; 
1  9  5  3) 
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where  r  and  0  are  sun  centered  polar  co-ordinates,  a  is  the 
aiming  parameter  (i.e.  the  radial  distance  of  the 
observation  point  from  the  center  of  the  solar  disk),  and  p 
is  the  real  part  of  the  complex  index  of  refraction.  Both  r 
and  a  are  measured  in  terms  of  the  photospheric  radius,  R  . 
The  optical  depth  is  related  to  the  absorption  coefficient, 
K,  hv  the  expression 
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where  As  is  the  path  length  t  li  rough  layer  m. 

m 

Tlie  absorption  coefficient,  k  ,  in  layer  ra,  is  a  function 

m 

of  the  plasma  frequency  ,  ^p,m  nntl  t*'e  c°llisiori  frequency, 

v  ,  in  this  layer,  and  the  radio  frequency,  f,  of  the 

c,m 

observations.  When  the  imaginary  component  of  the  refractive 
index  is  much  less  than  1  and  the  real  component,  p,  is 
close  to  1 ,  the  absorption  coefficient  can  be  approximated 
by  the  expression 


The  collision  frequency  is  the  sum  of  the  collision 

frequency  of  electrons  with  neutral  particles, v  and  the 
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collision  frequency  of  electrons  with  ions,v.  : 
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where  >^mis  the  density  of  neutral  particles,  'Je,mtlle 
electron  density,  which  is  assummed  to  be  essentially  equal 
to  the  ion  density,  and  T  is  the  plasma  temperature.  Mtcolet 
(1973)  gives  for  the  value  -  3.4  x  1 0“  1 0  ,  while 
Ginzburg  (1961)  gives  an  expression  for  in  the  range  of 
3  0-50  with  only  a  logarithmic  dependence  on  \'e  and  T.  In 
general  the  absorption  coefficient  also  depends  in  a  highly 
complex  way  on  the  magnetic  field,  but  for  the  frequencies 
of  interest  here  and  the  magnetic  field  strengths 
characteristic  of  coronal  holes,  the  magnetic  field  effects 
are  essentially  negligible,  and  can  therefore  be  neglected. 

Once  the  temperature  and  particle  density  profiles  are 
known,  the  ray  tracing  procedure  (equation  4)  can  be  used  to 
follow  the  path  of  the  ray  and  the  brightness  temperature  of 
the  emerging  ray  will  be  given  by  equation  2.  In  order  to 
account  for  different  changes  of  the  quiet  atmosphere  model 
In  different  height  ranges  inside  the  cornal  hole,  we  have 
divided  the  standard  atmospheric  model  into  five  regions: 
the  photosphere  (0  -  600  km),  the  inner  ch romosphere 
(600-2300  km),  the  outer  chromosphere  (2300  -  2465  km),  the 
transition  zone  (2465  -  2520  km),  and  the  corona  (2520  - 
700,000  km).  Since  our  radio  data  provide  no  information 
about  the  photosphere,  this  region  played  no  role  in  our 
analysis  and  therefore  was  left  unmodified.  The  upper  four 
regions  were  given  the  numerical  designations:  lower 
chromosphere  (0),  upper  chromosphere  (1),  transition  zone 
(2),  and  corona  (3). 

In  each  of  these  upper  four  regions,  the  electron 


density  of  each 
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density  can  he  modified  hy  multiplying  the 
layer  in  a  region  by  a  numerical  factor  (i»0,l,2,3). 

Similarly,  the  temperature  of  each  region  can  be  modified  by 
multiplying  the  temperatures  by  a  factor  8^  .  For  the 
modification  of  the  neutral  particle  density  we  use  a  single 
factor,  y,  for  all  the  regions  because  their  effects  are 
important  at  only  the  lower  regions  which  are  reached  only 
by  the  shorter  wavelengths.  A  different  form  of  modification 
is  achieved  by  changing  the  thickness  of  each  of  the  upper 
four  regions.  The  thicknesses  are  changed  by  multiplying 
them  with  the  factors  (i«  0  -  3),  which  in  effect  change 

the  gradients  of  all  quantities  in  each  region.  If  the 
modifications  in  two  adjacent  regions  are  different,  the 
corresponding  physical  quantities  are  interpolated  near  the 
hounJary  so  as  to  assure  a  smooth  transition  from  one  region 
to  the  next.  When  all  the  numerical  factors  mentioned  above 
are  equal  to  unity  we  simply  have  the  standard  solar  model. 

The  entire  procedure  of  modifying  the  standard 
atmospheric  model,  tracing  ray  paths  and  obtaining 
brightness  temperatures  at  several  wavelengths,  normalizing 
the  brightness  temperatures  to  the  corresponding  values  for 
the  unmodified  atmosphere,  and  finally  plotting  the  entire 
radio  spectrum  for  a  particular  modification  of  the  standard 
model  is  accomplished  by  a  single  computer  program.  This 
program  is  listed  in  Appendix  II  along  with  a  sample  input 
and  output. 

Although  it  is  possible  to  determine  the  effect  of 
changing  a  single  parameter,  in  general  these  parameters 


interact  with  each  other  in  non-linear  ways,  as  is  the  case 
between  electron  density  and  temperature.  \  particular 
example  is  the  relation  between  the  electron  density  and  the 
depth  of  the  transition  zone.  The  total  number  of  electrons 
in  the  transition  zone  does  not  appear  to  c  aan  .•  v.- r  . 
;■>  i  ",  n  i  f  i  c  >  n  t  1  v  is  it  expands  in  a  coronal  Vile.  "bus  when  'he 
transition  zone  expands  hy  a  factor  of  about  '>  or  (>  ,  is  was 
sn  ip’estod  by  XlIV  data  from  Sky  lab  (Bo  III  in,  et  al.,  1  9  7  '>  ; 
duber,  et  al.,  1974;  Withbroe,  1977),  the  electron  number 
density  decreases  to  about  20%  of  the  normal  value.  The 
effect  of  the  different  modifications  on  the  observed 
brightness  depends  also  on  the  aiming  parameter,  a,  since  at 
different  values  of  the  aiming  parameter  the  rays  follow 
different  paths  through  the  solar  atmosphere,  dear  the  limb, 
for  example,  at  longer  wavelengths  we  have  the  phenomenon  of 
limb  brightening,  which  tends  to  obscure  the  effects  of 
coronal  holes. 

It  must  he  mentioned  that  the  entire  modeling  procedure 
'..e  have  followed  assumes  a  spherically  symmetric  solar 
atmosphere  anJ  neglects  changes  toward  the  polar  regions, 
which  are  actually  present.  We  are  mostly  interested, 
however,  in  the  modifications  of  a  given  profile  of  the 
solar  atmosphere  inside  a  coronal  hole.  The  use  of  a 
somewhat  different  profile  In  the  polar  regions  would  simply 
produce  somewhat  different  modifications  more  appropriate 
for  this  model,  but  overall  the  changes  are  likely  to  be  in 


the 


s  a  me 


direction  as  in  the  case  of  the  standard  node  1 
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>  .  i  .  Mo d  v  1  i  in-  -  Results 

i’o  obtain  a  model  for  the  solar  atmosphere  inside  a 
coronal  hole,  we  begin  by  determining  the  effects  on  the 
spectrum  of  changing  each  parameter  of  the  standard  model 
separately.  These  effects,  of  course,  will  be  somewhat 
different  for  different  values  of  the  the  aiming  parameter, 
a.  if  changing  a  single  parameter,  such  as  or  T,  doesn't 
produce  t  lie  desired  results  we  try  cominations  of  the  two 
parameters  forming  a  g r i d  of  different  numerical  values  for 
at  and  ft.  The  fitting  of  the  coronal  hole  spectrum  proceeds 
from  the  longer  wavelength  end  because  this  corresponds  to 
the  outer  most  livers  and  is  not  affected  by  changes  in  the 
lower  regions.  We  then  proceed  to  the  shorter  wavelengths 
and  the  deeper  layers  of  the  atmosphere.  Once  a  reasonably 
good  fit  is  found  for  the  entire  spectrum,  the  parameters 
'•in  ho  further  modified  by  small  amounts  to  improve  the  fit. 

To  demonstrate  how  this  procedure  works  we  will  consider 
ia  detail  t  lie  process  as  it  applied  to  one  particular 
coronal  hole,  the  southern  coronal  hole  of  July  1973.  Figure 
*  >  shows  the  computed  brightness  temperature  ratios  and  the 
observed  ratios  for  the  southern  coronal  hole  of  July  17, 
19/3.  The  preliminary  estimates  for  the  parameters  are  given 
in  Table  13. 
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TABLE  10 

Preliminary  Parameter  Estimates 
Southern  Coronal  dole  Disk  Spectrum  -  July  17,  1978 
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In  Figures  49  -  56  we  show  the  procedure  used  to  refine 
the  original  estimates.  Each  graph  contains  4  curves  plus 
the  measured  spectrum.  In  each  case,  two  parameters  in  one 
region  (e.g.  ai  and  $1  )  have  been  modified  around  the 
initial  estimates.  From  this  set  of  graphs  the  effects  of 
varying  the  different  parameters  around  the  initial 
estimates  can  be  determined.  The  dominant  effect  of  changing 
parameters  in  the  chromosphere  (regions  0  and  1)  is  observed 
at  wavelengths  shorter  than  about  5  cm.  Varying  the 
parameters  in  the  transition  zone  (region  2)  primarily 
affects  the  wavelength  range  from  about  5  cm  to  50  cm,  while 
changing  the  corona  (region  3)  affects  primarily  the 
wavelengths  longer  than  50  cm  but  of  course  also  has  some 
affect  on  wavelengths  in  the  5  cm  to  50  cm  range. 

One  significant  effect  is  that  of  changing  the  thickness 
of  a  region  (parameter  c,  ).  Increasing  the  thickness  of  a 
region  while  keeping  the  density  and  temperature  unchanged 
results  in  an  enhanced  brightness  temperature  at  wavelengths 
which  emanate  primarily  from  the  modified  region.  However, 
if  at  the  same  time  the  density  is  decreased  so  as  to 
maintain  a  constant  number  of  particles  in  the  region,  the 
two  effects  tend  to  cancel  out,  and  only  a  small  change  in 
brightness  results. 
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After  observing  the  effects  produced  by  changing 
different  pairs  of  parameters,  improved  fits  to  the  observed 
spectrum  can  be  made  by  varying  all  the  parameters  over  a 
limited  range.  The  result  of  this  final  refinement  of  the 
parameters  is  shown  in  Figure  57  and  the  values  of  the 
parameters  are  given  in  Table  11. 
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Final  Parameter  Estimates 

Southern  Coronal  Hole  Disk  Spectrum  -  July  17,  1978 
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Our  best  results  for  the  coronal  hole  modeling  come  from 
the  July,  1978  observations.  We  have  used  these  data  to 
produce  models  of  both  the  southern  coronal  hole  and  the 
northern  coronal  hole  for  both  July  16  and  July  17.  In 
addition,  we  have  produced  separate  models  for  the  two 
different  aiming  parameters  for  the  southern  coronal  hole. 
These  results  are  summarized  in  Figures  57  -  62.  In  each 
figure  the  observed  spectrum  is  shown  as  open  boxes,  while 
the  theoretical  spectrum  is  shown  as  open  circles.  As  can  be 
seen,  the  models  fit  the  observations  quite  well.  The  values 
of  the  parameters  for  each  model  are  listed  in  Table  12.  We 
note  that  the  parameters  are  quite  consistent  both  for  the 
two  different  days  and  for  the  two  coronal  holes. 
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Southern  Coronal  Mole  Limb  Spectrum 
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Northern  Coronal  Hole 
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The  layers 

of  the  solar  atmosphere 

wh  ich 

c  on  t 

r  l  b  u  t  e 

the 

largest  fraction  of  the  radio 

flux  observed, 

vary 

f  ro  ;n 

wavelength  to 

wavelength.  Figures  63 

a  nd 

64 

show 

the 

contribution 

to  t  tie  flux  as  a 

funct ion 

of  height 

for 

different  wavelengths.  These  figures  are  derived  f ron  the 
model  of  the  southern  coronal  hole,  but  using  a  slightly 


different  modification  for  a  different  coronal  hole  does  not 
significantly  alter  the  contribution  profiles. 
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SOUTHERN  CORONAL  HOLE  LIMB  SPECTRUM  -  JULY  17 


SOUTHERN  CORONAL  HOLE  DISK  SPECTRUM  -  JULY  16 


NORTHERN  CORONAL  HOLE  SPECTRUM  -  JULY  16 
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AIMING  PARAMETER -0.70 


FIGURE  63 


CONTRIBUTION 
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AIMING  PARAMETER  *0.90 


FIGURE  64 


5.4 


Conclus  ions 


The  models  we  have  computed  for  coronal  holes  are  quite 
consistent  from  one  hole  to  another  and  from  one  viewing 
angle  (aiming  parameter)  to  another.  The  general  features 
are: 

1)  The  chromosphere  expands  in  thickness  by  about  a 

factor  of  2,  but  apparently  fills  in  with  material 
from  the  p  ho  t  os  ptie  r  e  ,  since  the  values  for  Nn  and 
both  decrease  by  only  about  2  0  %  .  The 

controversial  enhancement  near  2  cm  can  be 
obtained  through  a  significant  increase  of  the 
neutral  particle  density  in  this  region. 

2)  The  transition  zone  expands  in  thickness  by  about  a 

factor  of  5  to  6  and  the  electron  density  drops  by 
a  similar  factor. 

3)  The  corona  itself  does  not  appear  to  expand,  but  is 

depleted  in  particle  density. 

4)  The  temperatures  are  typically  reduced  by  about 
20%,  except  in  the  transition  zone,  where  the 
reduction  is  around  50%.  This  indicates  continued 
heating  of  these  regions  despite  their  expansion. 

The  coronal  hole  models  obtained,  together  with  other 
relevent  information,  such  as  solar  wind  speeds,  etc.,  can 
now  be  used  to  develop  a  quantitative  understanding  of  the 
physical  processes  operating  inside  coronal  holes  and  to 
bring  theory  and  observations  into  full  agreement.  This  is  a 
task  that  was  beyond  the  goals  of  this  project,  hut  we  hope 
to  be  able  to  undertake  it  in  the  near  future. 
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APPENDIX  I 

STANDARD  SOLAR  ATMOSPHERE  MODEL 
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08 

O.OOOOOe 

00 

87720.0 

1933300.0 

0,69391e 

08 

O.OOOOOe 

00 

77720.0 

1940400.0 

0. 75181e 

08 

O.OOOOOe 

00 

67720.0 

1947600.0 

0. 81603e 

08 

O.OOOOOe 

00 

57720.0 

1955000.0 

0, 88743e 

08 

O.OOOOOe 

00 

47720,0 

1962500,0 

0 . 9669  6e 

08 

O.OOOOOe 

00 

37720.0 

1970100,0 

0 , 15000e 

09 

O.OOOOOe 

00 

27720.0 

1959700,0 

0 , 15309e 

09 

O.OOOOOe 

00 

17720.0 

1692900.0 

0 ,17721e 

09 

O.OOOOOe 

00 

16720.0 

1659900.0 

0,18074e 

09 

O.OOOOOe 

00 

15720.0 

1625200.0 

0 , 18460e 

09 

O.OOOOOe 

00 

14720.0 

1588600.0 

0.18385e 

09 

O.OOOOOe 

00 

13720,0 

1549700.0 

0. 19358e 

09 

O.OOOOOe 

00 

12720.0 

1508400,0 

0 , 19889e 

09 

O.OOOOOe 

00 

11720.0 

1464100.0 

0, 20491e 

09 

O.OOOOOe 

00 

10720.0 

1416200,0 

0 , 21183e 

09 

O.OOOOOe 

00 

9720.0 

1364000.0 

0, 21994e 

09 

O.OOOOOe 

00 

8720.0 

1306400,0 

0. 22964e 

09 

O.OOOOOe 

00 

7720.0 

1241800,0 

0, 24158e 

09 

O.OOOOOe 

00 

6720.0 

1167700,0 

0 . 2  569 le 

09 

O.OOOOOe 

00 

5720,0 

1080000,0 

0, 27779e 

09 

O.OOOOOe 

00 

4720.0 

970070,0 

0 , 309  26e 

09 

O.OOOOOe 

00 

3720,0 

816610,0 

0, 36737e 

09 

O.OOOOOe 

00 

2720.0 

503800,0 

0, 59 548 e 

09 

O.OOOOOe 

00 

2620.0 

428270.0 

0 , 70049e 

09 

O.OOOOOe 

00 

2520.0 

286250.0 

0 , 1048  Oe 

10 

O.OOOOOe 

00 

2510.0 

257990,0 

0 , 116  28e 

10 

O.OOOOOe 

00 

2500,0 

219070,0 

0 , 1369  4e 

10 

O.OOOOOe 

00 

2495.0 

190820,0 

0,15722e 

10 

O.OOOOOe 

00 

2490.0 

145250,0 

0. 20654e 

10 

O.OOOOOe 

00 

2489.0 

130290,0 

0, 23025e 

10 

O.OOOOOe 

00 

2488.0 

109290.0 

0, 27450e 

10 

O.OOOOOe 

00 

2487.0 

65000,0 

0 , 69  420e 

10 

0.99942e 

06 

2485.0 

45000.0 

0 ,97000e 

10 

0, 59965e 

07 

2484.0 

34000,0 

0, 12390e 

11 

0, 3999  7e 

08 

2483.0 

27000,0 

0 , 149  20e 

11 

0, 21000e 

09 

2482.0 

24000,0 

0, 15770e 

11 

0. 52000e 

09 

2481,0 

22500,0 

0 . 15790e 

11 

0, 83000e 

09 
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HT(km) 

T(°K) 

N(e1 /cm3 ) 

Nn(part/cm3 ) 

2480.0 

22000.0 

Q,15820e 

11 

0 , 96GQ0e 

09 

2465.0 

21000,0 

0 , 16420e 

11 

0 . 10  2  00e 

10 

2445.0 

20600,0 

0, 17090e 

11 

0 ,9  7000e 

09 

2415.0 

20400,0 

0. 17830e 

11 

0 ,94000e 

09 

2375.0 

20300.0 

0 , 18660e 

11 

0,9  70Q0e 

09 

2315.0 

20200.0 

0, 19670e 

11 

0 . 11600e 

10 

2305.0 

19800.0 

0 , 20110e 

11 

0, 13900e 

10 

2302.0 

16000,0 

0,22670e 

11 

0 , 33700e 

10 

2298.0 

13000,0 

0 , 2  50  20e 

11 

0 , 63600e 

10 

2295.0 

8800,0 

0 , 2  899  Oe 

11 

0. 13930e 

11 

2282.0 

8550,0 

0 , 28  360e 

11 

0. 17760e 

11 

2268.0 

8300.0 

0, 2  7600e 

11 

0 , 21880e 

11 

2235.0 

7800,0 

0, 26260e 

11 

0, 32230e 

11 

2170.0 

7300,0 

0, 24590e 

11 

0. 52880e 

11 

2095.0 

6940,0 

0 , 229  70e 

11 

0. 84160e 

11 

2030.0 

6720,0 

0 , 218  40e 

11 

0 . 12207e 

12 

1945.0 

6620.0 

0 , 22210c 

11 

0, 18554e 

12 

1810.0 

6560.0 

0, 25280e 

11 

0. 35890e 

12 

1690.0 

6500.0 

0, 29320e 

11 

0. 66508e 

12 

1560.0 

6440.0 

0 , 3  748  Oe 

11 

0 , 1352  6e 

13 

1430.0 

6390,0 

0, 51500e 

11 

0. 28774e 

13 

1280.0 

6330.0 

0, 79380e 

11 

0 , 72643e 

13 

1170.0 

6230.0 

0 ,l0550e 

12 

0 . 148  50e 

14 

1030.0 

5930,0 

0. 12690e 

12 

0. 39885e 

14 

980.0 

5760.0 

0 , 12 180e 

12 

0. 58694e 

14 

930.0 

5560,0 

0,11060c 

12 

0 , 88 170e 

14 

880.0 

5360.0 

0.97280e 

11 

0.13413e 

15 

830.0 

5150.0 

0 , 85560e 

11 

0.20715e 

15 

780.0 

4890,0 

0 , 77910e 

11 

0. 32997e 

15 

730.0 

4600,0 

0.84240e 

11 

0, 54269e 

15 

680,0 

4350,0 

0, 11220e 

12 

0 , 91350e 

15 

630.0 

4170.0 

0, 16150e 

12 

0 . 15560e 

16 

605.0 

4150,0 

0, 19  610e 

12 

0 , 20060e 

16 

580. C 

4150,0 

0, 24100e 

12 

0 , 25760e 

16 

530.0 

4200,0 

0, 37540e 

12 

0 . 418  70e 

16 

480.0 

4350.0 

0, 61510e 

12 

0,65730e 

16 

430,0 

4500,0 

0 , 99040e 

12 

0 , 10160e 

17 

380.0 

4640.0 

0, 15620e 

13 

0, 15540e 

17 

330.0 

4780.0 

0, 24310e 

13 

0 . 23460e 

17 

280.0 

4930,0 

0, 37610e 

13 

0 . 349  40e 

17 

230.0 

5100,0 

0, 58460e 

13 

0. 51170e 

17 

180.0 

5370,0 

0,96450e 

13 

0. 72439e 

17 

130.0 

5750.0 

0 , 18980e 

14 

0,98574e 

17 

100,0 

6070.0 

0, 348  50e 

14 

0 , 11548e 

18 

80.0 

6300,0 

0 , 55110e 

14 

0, 12736e 

18 

65.0 

6600,0 

0 , 95490e 

14 

0, 1340  3e 

18 

50.0 

7100,0 

0, 22090e 

15 

0. 13641e 

18 

30.0 

8000,0 

0, 80400e 

15 

0. 13433e 

18 

0.0 

9000,0 

0 ,26010e 

16 

0. 13495e 

18 
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APPKNDIX  II 
RAY  TRACING  PROGRAM 


/ 1  oad 
/  pa  rm 
/job 
c 

cccc 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


1  4 
7  1 


c 

c 

cccc 

c 

c 

c 

c 


f  o  r  t  g  1 , n  s  egs  =6 
deck 
nogo 


program 
d ime  ns  ion 


"solmfw" 


dh(200) , h( 200 ) , t ( 200 ) , c  r ( 200 ) ,dtau(200)  ,  c  rm ( 2  l  )  , 
ar(200),br(200),dr(200),bt(21),btq(21),spec(21) 
r,rmml,s,com(40),rcrm(2l),rra,ocr(51),rl,r2,r3,rlpl, 
r2pl  ,r3pl  , p , q , r mp 1  ,ppl  ,pp2,rnl  ,  r  t ,  r  3n  1  , r 4 , r4n 1  ,r0, 
rOpl , rlm3 

l(2l),nu(21)  , nn(200) ,mu2( 200) ,  1  nu  ,  ne(  200)  ,  htel  (  2  1  )  , 
hcrl(21),hcr2(21),nue(200),nuep(200),nuen(200), 
tq(200) ,ncq( 200) , ro2( 200) , hr(21 ) , tauref(21 ) , the(200) , 
the inf ( 2 l ) ,er(200) ,xr(200) , nnq(200) , hq( 200) 
logical  qflag/. false. / , a  f 1 g / . false./ 
integer*4  gtitl(15) 


integer 


real 


r m» 1 8  1 

rm  «  no.  of  levels  in  the  solar  atmosphere, 
r  mp  1  -  rm+  1 
rmml -rm- 1 
xc-2 . 9  9  7  9e  1  0 

xc-  speed  of  light  (cm/s), 
r sun-6 . 9  59  9e  5 

rsun  *  radius  of  the  sun  (km). 
pi-3. 14] 5927 
pi  -  that's  right,  pi. 
crd-180.0/pi 

crd  -  conversion  of  radians  to  degrees. 
x5- 1 .0e5 

x5  *  number  of  cm/km. 
xk  f c- 1 . 1 e-6 

xkfc  »  constant  for  computing  conductive  flux. 
xh3-3 . 13260e9 

xb3  -  constant  for  computing  the  plasma  frequency. 

do  14  1-1,51 
oc  r ( i ) -0 
do  71  i « 4 , 3 6 , 4 
ocr( i  +  1 )  —  1  /  4 
ocr ( 4 1 ) - 1 
x9- 1 .0e-9 


input  wavelengths,  aimpar,  spectrum,  and  region  limits  cccccccccc 
read  ( 5  ,  * )  im 


lm  is  the  number  of  wavelengths  that  are 


to  be  processed 


i 

1 


i 


do  16  i-1  ,  ira 


1  28 


16  read  (5,617)  l(i) 

617  format  (3fl0.0) 

c  l(i)  “  wavelength  (cm). 

c  im  =  number  of  wavelengths  to  process. 

HOI  read  (3,602,end*998)  almpar,spec 

602  format  ( f 7 . 0 , 2 1 f 3 . 2 ) 

read  (5,851)gtitl 
8  5  1  f  o  rina  t  ( 1  5aA  ) 
a  f 1 g“ . false. 

c  spec(i)  “  coronal  hole  spectrum. 

d2»aimpar*aimpar 

read  (5,611)  r(),rl,r2,r3,lp,alfaO,alfal,alfa2,alfa3,betaO, 
l  be t a 1 , be t a 2 , be t a3 , gamma 

611  f  o  rma t(515,9f6.0) 

c 

c 

i  f ( q  f 1 a  g ) g  o  to  30  2 

read  (5,601)  (  coir  (  l  )  ,  l  *  1  ,  AO  ) 

qflag*.true. 

601  format  (A0a2) 

c 

do  1  r  » 1 , r  m 

read  (  5  ,  6 00  ,  e nd  =  9 93  )  hq( r ) , tq( r ) , ncq ( r ) , nnq(r) 

600  format  ( f 8 . 0 , f 1 0 . 0 , e 1 2 . 0  ,  e 1 1  .  0 ) 

c  liq(r)  *>  altitude  of  layer  r  (km), 

c  tq(r)  -  temperature  in  layer  r  (k). 

c  ncq(r)»  <1  electrons  per  cm**  3 

c  nnq(r)»  !t  neutral  atoms  per  cm**  3 . 

1  continue 

c 

302  a  If  a0»  1  .0 

alfal-1  .0 
a  1 f  a2“  1  .0 
a  1 f  a3  *  1 .0 
b  e  t  a 0 *  1  .0 
betal  =*  l  .0 
b  e  t  a 2  *  1 .0 
b  e  t  a  3  *  1  .0 
ze  ta0»  1  .0 
ze  ta  1  ■  1  .0 
ze  ta2*  1  .0 
ze  ta3»  1  .0 
g  a  m  m  a  « 1  .0 
lcount“0 

the(rmpl)»crd*arsln(airapar*rsun/(h(rmpl)+rsun)) 

c 

cccc  modify  the  quiet  sun  model  c c ccc ccccccccc c ccc ccc ccc cc cccc cc cccc c c 
c 

7  r0pl-r0+l 
r  1  p  1  ■  r  1  +•  1 
r2  pl»  r  2+ 1 
r  3  p  1  *  r 3+ 1 
r3m 1 ■ r  3-  1 
r A- r 3- 2 


-  129  - 


r  Am  l  <*  r4-  1 
r  1  m  3*  r  1  -  3 
do  9  r*l  ,rO 
t ( r )»tq( r  ) 

9  ne(r)»neq(r) 
c 

t  p*  t  q  (  rO  ) 

do  6  r  »  r  0  p  1  ,  r  1 

t  (  r)*t  p+beta()*(  tq(r)-tp) 

6  ne(r)=neq(r)*alfaO 
c 

t  p«*tq  (  rO) 

do  8  r=rl pi , r2 

t(r)=tp+betal*(tq(r)-tp) 

if  ( t ( r ) . gt . t q ( r 2 ) )  t(r)»tq(r2) 

3  ne(r)»neq(r)*alfal 
c 

t p-t ( r2  ) 

Jo  19  r=*  r  3  ,  rm 
t(r)=tp+beta3*(tq(r)-tp) 

19  ne( r ) = neq ( r )*a 1 f a3 
c(rmpl )>t(rm) 
c 

t  p*  t ( r 2  ) 

do  17  r  *  r  2  p  l  ,  r  3ral 
t(r)»tp+beta2*(tq(r)-tp) 

17  ne ( r )= neq ( r ) *a If a2 

c 
c 

c  recompute  the  altitude  of  each  level 

c  this  also  changes  the  gradients  in  each  region 

c 

do  900  r = 1 , rO 

900  ti(r)*hq(r) 
hl»hq ( rO) 

do  9QL  r=r0pl,rl 

h(r)»hl  +  zetaO*(hq(r)-hq(r-l)  ) 

901  hl-h(r) 

do  902  r»rl pi , r2 
h(r)*hl+zetal*(hq(r)~hq(r-l)) 

902  hl-h(r) 

do  903  r*r2pl,r3ml 
h(r)»hl+zeta2*(hq(r)-hq(r-l)) 

903  hl-h(r) 

do  904  r^rS.rm 

h(r)»hl+zeta3*(hq(r)-hq( r-1  )  ) 

904  hi »h ( r  ) 

do  905  r*l  ,  ria 
x*(rsun+h(r))/rsun 
ro2 ( r ) »x*x 

905  cont  tnue 
c 

do  2  r»2 , rmm 1 
dh(r)-(h(r+l)-h(r-l))*0.5 
dh(l)-h(2)-h(l) 


2 


1  30 


c 


2  3 
c 
c 


2  9 
c 
c 


2  (> 
c 


c 


c 

CCCf 

c 

c 

c 

c 


d  h  (  rm  )  *  h  (  r  m  )  —  h  (  r  mm  1  ) 
dh(r)=  shell  thickness  (km), 
li (  r in p  1  )»h(rm)+0.5*dh(rm) 

natch  the  electron  density  at  r4**>r3. 

s  3  =  (  n  e  (  r  3  p  1  )-ne(  r3)  )/(ti(r3pl  )  -  h  (  r  3  )  ) 

do  23  r “ r  2  p 1 , r  3 

nep  =  ne(r3)+s3*(h(r)-h(r3)  ) 

if  ( ne ( r ) . ge  .  nep )  go  to  23 

no( r)=nep 

contl nue 

match  the  temperature  at  r4==>r3. 
s4  =  ( t(r3)-t(r4) ) / ( h ( r 3 ) -h ( r 4 )  ) 
do  20  r=r4,r3 

t( r)  =  t(r4)+s4*(h(r)-h(r4)  ) 

natch  the  electron  density  at  r * r 1 m 3==> r 1 + 1 . 
sl=(ne(rlpl)-ne(rlm3))/(h(rlpl)-h(rlm3)) 
do  26  r*rlm3,rlpl 
ne(r)*ne(rln»3)+sl*(h(r)-h(rlni3)) 


do  11  r  =  1  ,  r m 
gam* g  a  m  m  a 

if  ( gam . eq . 1 . 00  )  go  to  11 
if  ( n nq ( r  )  .  le  .  0 . 0  )  go  to  11 
a  If  ae  =  ne(r)/neq(r) 
hetae*tq(r)/t(r) 
betae  is  inverted! 

ga:nmax*betae+(  2 .0*neq(r)/nnq(  r  )  )*(betae-alfae) 
if  ( gamma . gt . gammax)  gam*gaminax 
nn(r)=nnq(r)*gam 


electron  collision  frequencies 


ICCCCOCCC 


nlcolet  ,  in.  (1  0  53),  "the  collision  frequency  of  electrons  in  the 
Ionosphere",  journal  of  atmospheric  and  terrestrial  physics, 
v  o 1 .  3,  pp.  200-2  1  1  . 


c  const,  nts  for  collisions  with  neutrals. 

21  xkn  =  2 . 6  0  2  2  e  6 

s  1  g  =  2 . 0  e  -  8 

c  sig  *  collision  distance  for  electrons  with  neutrals  (cm). 

sig2*sig*sig 

c  constants  for  collisions  with  positive  ions. 

x32-(-3. 0/2.0) 
x2 3- (-2. 0/3.0) 
x  b I  *  I  .81668 
xk4-5 . 72965e6 
c 

do  15  r  *  1  ,  r  m 

nuen(r)«xkn*nn(r)*sig2*sqrt(t(r)) 

c  nuen(r)  *  collision  frequency  between  electrons  and  neutrals. 


al2*alog( l . 0+(xk4*t(r)*t(r)*(2.0*ne(r))**x23)  ) 


1  3  1 


c 

c 

1  5 

c 

c 

c  c  c  c 
c 


003 


c 

3  1  l 


0  04 

613 

22 

c 

c  c  c  c 
c 

c  c  c  c 


c 


5 

c 

cc 

cccc 
c  c 
c 


c 


nuep( r ) »x  b 1 *a  1 2*ne( r  )  *  ( t( r)**x32) 

nuep(r)  ■=  collision  frequency  between  electrons  and  protrons. 
nue( r ) -nuep ( r)+nuen(r) 

nue(r)  »  electron  collision  frequency. 

i f ( 1 p . eq . -  1 ) g o  to  311 

write(ll,616)(com(j),j»l,40),r0,alfa0,beta0,zeta0,rl,alfal, 

1  betal,zetal,r2,alfa2,beta2,zeta2,r3,alfa3,beta3,zeta3. 


2  ga  mma , a i mpa  r 
write  (11,603) 
f  o  rma  t  ( 

1  0  5  h  r 

h 

dh 

t 

ne 

1  n  n  nue 

nu  e  p 

nue  n 

ratio 

f 

p*p/f n 

2 , 1 0  2  h  (-)  (km)  (km)  (k)  ( c  m*  *- 3 )  (cm**-3)  (hz) 

3  ( h  7. )  (Hz)  c/s*cm*cm  /) 

do  22  r ■  2  ,  r in 
ratio=nuen( r)/nuep( r) 

fcond“xkfc*t(r)*t(r)*sqrt(t(r))*(t(r+l)-t(r-l))/(h(r+l)-h(r-l))‘ 
f condd*f  cond 

if  (fcondd.eq.0.00)  fcondd*1.00 

p2of»ne(r)*t(r)*ne(r)*t(r)/fcondd 

alfae=*ne(r)/neq(r) 

be  t ae* t ( r ) / 1  q ( r ) 

i f ( 1 p . eq . -  1 ) go  to  22 

write  (  1  1  ,604)  r ,  h ( r ) , dh ( r ) , t ( r ) , ne ( r ) , nn ( r ) , nue ( r ) , nue p ( r )  , 

1  nuen(r),ratio,fcond,p2of,alfae,betae 

format  (  i  6  ,  £8 . 0  ,  f  8 . 1  ,  f  9 . 0 , 1  |>5e  1  0 . 3 , 5e  9 . 2  ) 

if  ( r . eq . r 1 . or . r . eq . r 2 . o r . r . eq . r 3  . o r . r . eq . rO)  write  (11,613) 

format  (lx,13(10h . )) 

cont  inue 

optical  depth  and  contribution  function  routines  cccccccccccccccc 

initialize  variables, 
do  5  i”l  ,  im 
nu(i)*xc/l(l) 

nu(i)  ”  frequency  corresponding  to  l(i)  (hz). 

bt(i)-0.0 

htel ( i )«0.0 

her  1 ( i )  =  0.0 

hcr2( i  )*0.0 

rcrm( i  )*l 

crm(i)**0.0 


wavelength  loop. 

do  999  i-1  ,  im 
lnu-a log( nu ( i  )  ) 

if  (l(i).ge. 0.030)  go  to  13 


**•-**=. 


(>  n 


write  (6,607)  (com(j),j*l,40),l(i) 

write  (6,613) 

format  (5x,  57hthe  r e y 1 e i gh- j e a  ns  approximation  is  used  in  this 
1 rograrn . / 5x ,  72hthis  makes  computations  at  wavelengths  less  than  0 
2050  cm  inappropriate.) 
go  to  999 
c 

cccc  compute  dtau  for  each  layer  and  locate  the  turning-point  level. 

13  ltl=l(i)*l(i) 

0171*2. 0*pi*nu(  i  ) 

om2*om*om 

mu2 ( rnp 1 ) *  1  .00 

p  =  0 

c 

do  10  j  = 1 , r m 

r  =  rmp 1  - j 

omo2  =  xb3*ne(  r) 

yp*om2+nue(r)*nue(r) 

xm f  *  1 .00-(omo2/yp) 

xnf  =  nne( r ) *  o  m  o 2 / (om*yp) 

ypp*(xnf*xnf)/(xmf*xmf) 

chi2=0.25*xnf*ypp*( 1 .0-0.25*ypp) 

mu2(r)“xmf+chl2 

if  (mu2(r).le.0.00)  go  to  60 
c  if  (mu**2.1e.0.00)  the  wave  is  reflected, 

z  i*  1  . O0-(d2/(mu2(r)*ro2(r))) 
if  (zi.le.0.00)  go  to  60 

c  if  (zi.le.0.00)  the  turning-point  has  been  reached. 

dtau(r)*2 .O*sqrt(chi2)*om*x5*dh(r)/ (xc*sqrt(zi)  ) 
if  (dtau(r)-l75.0)  25,24,24 


24 

a  r ( r )  =  1 .00 

go  to  10 

2  5 

a  r ( r )  =  1 .00- 

exp(-dtau(r)) 

10 

con  t  i  nu  e 

go  to  61 

6  0 

p=r 

c 

p  =  pointer 

to  the  level  which  reflects  the  radio  wave. 

c 

C  C  C  r 

(’  c).  1  p  II  l  •!  ;)  t 

icnl  i opt hs  of  direct  and  reflected  rays. 

6  1 

p  P 1  =  !>  H 

d  r( rmpl ) *0.00 
dtau(rmpl)*0.00 
j  t  *  rmp 1 - p 
do  30  j  =  l  ,  jt 
r  =  rin  pi  -  j 
r  p  1  *  r+ 1 

30  dr(r)*dr(rpl)+dtau(rpl) 
c  dr(r)  *  optical  depth  of  level  r  (direct  ray). 

62  pp2*p+2 

do  27  r*pp2,rm 
e  r ( r )*dr( p) 
rml*r-l 

do  27  q*ppl , rml 
27  er(r)*er(r)+dtau(q) 

er(r)  -  optical  depth  of  level  r  (reflected  ray). 


c 


-  133- 

er(ppl )  -  d  r  (  p  ) 
t.iuref  (  l  )  =  !  r  (  p) 

c  tauref(i)  =  optical  depth  of  turning-point, 

c 

cccc  compute  contribution  from  each  level, 
do  32  r  =  p  p  1  ,  r  m 
b  r  (  r  )  =  0 . 0 

if  (dr(r)-l75.0)  31, 2d, 2 .3 

31  br(r)*hr(r)rcxp(-dr( r)  ) 

2  3  if  ( e r ( r ) -  1 7  3 . 0  )  3  6,32  ,  32 

36  br(r)=br(r)+exp(-er(r)) 

32  x r ( r )  =  t ( r ) *a r ( r ) *b r ( r  ) 

c  xr(r)  =  contribution  from  level  r  (k). 

c 

cccc  find  the  brie,  t.  temp.,  t  he  peak  of  t  lie  cont.  fn.  and  its  location, 
do  4  0  r  =  p  p 1  , r m 
cr(r)=xr(r)/dh(r) 

c  cr(r)  =  contribution  function  (k/km  at  level  r). 

If  (cra(  i  )-cr(r))  33,40,41) 

33crm(i)=cr(r) 

c  crm(i)  =  maximum  value  of  cr(r). 

r  c  rm  (  i  )  =  r 

c  rcrn(i)  =  location  of  the  peak  value  of  cr(r). 

4  0  b  t ( 1 )  =  b  t ( i )  +  x  r  (  r ) 

c  bt(l)  =  brightness  temperature  at  wave  le  light  1  (  i  )  (k). 

c 

cccc  normalize  the  contribution  and  the  contribution  function, 
c 

if  (cr m (i).eq. 0.00)  go  to  72 
x  a  *  1  . 0  /  c  r  m  (  I  ) 
do  41  r=ppl , rm 
41  cr(r)=>cr(r)*xa 

c  cr(r)  is  now  tbe  normalized  contribution  function. 

c 

72  xa-0.0 

do  6  9  r ” p p l , r ra 
if  (xa-xr(r))  63,69,69 
0 i  xa*xr(r) 

69  conti nue 
c 

if  (xa.eq.0.00)  go  to  73 
xa  » 1  . 0 / xa 
do  67  r-ppl , rm 
67  xr(r)*=xr(r)*xa 
c 

cccc  define  variables  for  output. 

7  3  ,ir(rnpl  )*0.D'l 
b  r ( rm  p 1 ) -  1 . 00 
c  r ( rmp l ) »0 . 00 
er(rmpl  )  *»  0 . 0  0 
xr(  rtnpl  )»0 .00 
d  t  au  ( |i  )  »0 . 0 
ar(p)-O.Oi) 
br ( p  )  *0 . 00 

c  r  (  p  ) «  0 . 0  0 


1  3  4 


c 

c 

cccc 


/ 1) 
c 

c 

cccc 
4  9 


4  6 

c 

74 

c 

4  7 

4  a 

c 

c 

>') 

5  1 
c 

52 

5  a 

5  4 

5  5 

c 

c 

cccc 


1  t  p  )  = ;) .  0  0 

x  r  (  p  )  =')  .00 
L  he ( p  )  =0 . 0 
h  r  (  i  )  =  h(  p  ) 

Ur(i)  =  licit*  lit  of  reflection  level. 

compute  theta  at  each  layer. 

y=*aimpar*crd/r8un 

r  t =  rm  hp  p 1 

do  70  j  =  p  p  1  ,  r  m 

r=  r t - j 

x=*y*dh(r)/sqrt(ro2(r)*(nui2(r)*ro2(r)-d2)) 
t  he  (  r  )  « t  lie  (  r-H  )+x 
cont  i n  u  e 

tho(rnpl)  was  defined  above. 

1  he  i  n  f  (  i  )  *  t  lie  (  p  p  1  ) 

determine  the  altitudes  of  the  half  power  points  of  cr(r). 
if  (  c  r  n  (  i  )  .  e  >| .  0 . 0  0  )  go  to  74 
j  *  r  c  r  in  (  i  ) 

.1  m  1  =  .)-' 

do  4  5  r  =* n  ,  j m  I 

«•*  i  +  l>-  I"  r 

if  ( c r ( s ) - 0 . 5 )  4  0,47,46 

cont i nu  e 

her 2  (  i  )  =  (-)99 .0) 

ho r2 ( i )  =  -999 .0  means  that  the  model  does  not  go  deep  enough. 

■go  to  50 

-i  c  r  2  (  i  >  =  (-32.3.0) 
hcrl  (  l  )  =  (-3fM.O) 

h  o  r  •'  (  i  )  =  -3 ,3.3.0  means  there  is  no  contribution  at  this  wave  Ion  t  n . 

o  to  5  6 
lie  r2  (  i  )  =  h(  s  ) 

'•O  to  5  0 
r  =  s  1 1 

h  c  r  2 (  i  )  =  ((0.5-cr(r))*(h(s)-h(r))/(cr(s)-cr(r)))  +  li(r) 

Iicr2(i)  =  lower  half-power  point  of  the  cont.  fn. 

if  (j-rmml)  52,51,51 
he  r  1  (  i  )  =  (-'),' 9  .  0  ) 

lie  r  1  (  i  )  =  -9  9  9 . 0  means  that  the  model  does  not  p,o  high  enough. 

>*  o  to  5  6 
Jo  55  r  =  j , r  m  p 1 
s  ---  r 

if  ( c  r ( s  )  - 0 . 5  )  5  5,54,53 

cont  In u v 

licrl  (  i  )  =  (-  999.0) 
go  to  56 
licrl(  i )=h ( s  ) 
go  to  56 
r  =  s-  1 

hcrl  (  i)  =  ((0.5-cr(r))*(h(s)-h(r))/(cr(s)-cr(r)))  +  li(r) 
hcrl(i)  =  upper  half-power  point  of  the  cont.  fn. 

output  contribution  function  for  wavelength  l(i)  cccccccccccccccec 


56  if  (ln.no. 1)  go  to  6  4 

57  wri  te(  1  I  ,616)(com(j),j=l  ,  4  0  )  ,  r  0  ,  a  1  f  a  (1  ,  b  o  t  a  0  ,  z  e  t  a  0  ,  rl  ,  a  1  f  a  1  , 

1  b  e  t  a  I  ,  z  e  t  a  1  ,  r  2  ,  a  1  f  a  2  ,  b  e  t  a  2  ,  z  e  t  a  2  ,  r  3  ,  a  1  f  a  3  ,  b  e  t  a  3  ,  z  e  t  a  5  , 

2  gain  tin  ,  a  imp  a  r 

urite  (  1  1  ,  607  )  1  (  1  ) 

607  format  (5x/16h  wavelength  =,19.3,56  ( c m ) , 6 4 x , 2 5 he on t r i hu t 1  on  f 
lrom  level  r  /  / 

1  13  )li  r  6  rrn  mu*  *2  dtau  ar  dr  cr 

2  br  theta  0.0  0.2  0.4  0.6  i.H  1. 

30  / 

4  1301,  (-)  (km)  (k/km)  (---)  (--)  (--)  (--)  (-- 

5)  (--)  (dee,)  ~ * ~ - 

6  /) 
c 

do  34  r=p,rmpl 
kmax»xr(r)*4 0.0+1 .001 
if  ( k  m  a  x . g t . 4  1  )  k  m  a  x  =  4  1 

write  (  1  l  ,603  )  r  ,h(r  ) ,  cr(  r)  ,i.iu2(  r  )  ,dt  au(r)  ,ar(  r  )  ,d  r(  r)  ,  e  r  (  r  )  , 

1  br(r),the(r),(ocr(k),k=l,kmax) 

60J  format  (  i  5  ,  f  3  . 0  ,  f  7  .  4  ,  f  7  . 4 , 5  e  1  0  . 3  ,  f  7  .  2  ,  2  x  ,  4  1  i  1  ) 

if  (r.eq.rl.nr.r.eq.r2.or.r.eq.r3.or.r.eq.r0)  write  (11,615) 

34  continue 
c 

cccc  find  the  height  at  which  tau  equals  unity. 

44  rt* rn- p 

do  39  r= 1  ,  r t 
s=*rni-r 

if  (dr(s)-l .0)  39,42,43 

39  continue 

htel(i)=(-999.0) 

c  htel(i)  =  -999.0  means  tau*l  not  reached  above  turning-point, 

go  to  999 

42  htel(i)=h(s) 
go  to  999 

43  r =s+ 1 

Intel  (  i  )  =  (  (  1 .0-dr(r)  )*(h(s)-h(r) )  /  (dr(s)-dr(  r)  ))+h(  r) 
c  htel(i)  =  In  eight  at  tau  =  l  at  wavelength  1(1)  (kin). 

999  continue 
c 

cccc  wavelength  dependent  parameters  ccccccccccccccccccccccccccccccccc 
c 

if(ln,eq.—  l);}o  to  312 

wr  i  t  e  (  1  1  ,616)(com(j),j“l  ,  40  )  ,  r  0  ,  a  1  f  aO  ,  be  t  aO  ,  ze  t  aO  ,  r  1  ,  a  1  f  a  1  , 

1  bet a l  ,  z  e  t  a  1  ,r2,alfa2,beta2,zeta2,r3,alfa3,beta3,zeta3, 

2  gamma , a  imp a r 
312  do  997  i  “  1  ,  i r.i 

nu(i)=»nu(i)*x9 
s»rc  rm(  i) 

Inu »alog 10(1(1)) 
if(lp.eq.-l)go  to  310 

997  write  (11,610)  l(i)tnu(l),bt(l),lnu,h(s),htel(l),crm(i),hcr2(i), 
l  licrl(l),hr(l),tauref(i),thelnf(i) 

format  (  5  x  ,  l  4  hwa  ve  l  e  ng  t  6  =  ,  0  p  f  1  0 . 3  ,  2  6 1,  cm  *  frequency 

1»,  f  9. 3, 2  5b  ghz  *  brlgt.  teinp.  =,  fl().(),7h  k  */ 


6  10 


-  Ilf)  - 


c 

;<  l  o 

c 


f,  \ ') 


c 


c 

20 


3  0  '3 


c 

c 

c 

c 

c 


3 

4 


c 


3x,l4hlog  u.ivlngt  .  =»  ,  fl  0.3,  2  fill  —  *  nc.  cfn.  pk. 


3-,  fO.l 

,  2  3  li 

It  m  *  h  t  . 

t  a  u  «  1.0  = 

,  f  1  0  .  1  , 

7  h 

k,n 

*  / 

4 

3  x , 

1 4  h  cfn. 

pit 

.  =  ,  1  p  e  1  0 . 3 

,  2  f>  h 

k  /  k  m 

* 

hts 

.  cfn.  =  0,3 

3= ,0pf 9  .  1 

,4b 

=  =  >  ,  f  9 . 1  ,  3  b 

km 

,  2  3  x ,  1 h*  / 

() 

3  x  , 

l  4  li  re  f  1  e  c  t . 

lit 

.  =  ,  f  1  0 . 1 

,  2bb 

km 

* 

t  a  u 

at  reflect. 

7= , 1  pc  9 . 2 

,2  3b 

_  * 

t  h . 

reflect.  = 

,  Op f 10.3, 

7  h 

tl  G  g  • 

*/) 

if  ( l count . eq . 0 )  go  to  3 

w  r  i  t  e  (  f>  ,  0  1  f>  )  (  c  mi  (  j  )  ,  j  =  l  ,  4  0  )  ,  r  0  ,  n  1  f  aO  ,  be  t  a  0  ,  ze  t  a  0  ,  r  1  ,  a  1 1  a  1  , 

1  be  t  a  1  ,  ze  t  a  1  ,  r  2  ,  a  1  f  a 2  ,  be  t  a 2  ,  ze  t  a 2  ,  r  3  ,  a  l  f  a  3  ,  be  t  a 3  ,  z e  t  a  3  , 

2  g  a  :n  in  a  ,  a  i  m  p  a  r 

format  (  1  li  1  /  A  x  ,  4  Da  2  /  /  3x  ,  f»h  r  0  =  ,  i  3  ,  1  1  h  ,  a  1  f  a  0  =  ,  f  3  .  2  ,  1  1  li  ,  be 
It  an  =  ,  f  5  .  2  ,  1  1  h  ,  zetaO  ”  ,  f  5  .  2  /  3  x  ,  6  h  rl  =  ,  i  3  , 

1  1  1  li  ,  a  1  f  a  1  =  , 

'»  f3.2,llh  ,  beta!  =  ,f3.2,lili  ,  zetal  =  ,f3.2 

1/  3x,i>li  r>  =,i3,llh  ,  alfa2  =,f5.2,llb  ,  beta2  =,fi.2, 

.lib  ,  zeta2  =,f5.2/  3x,6h  r3  * , i 3 , 1 1 h  ,  alfa3  =,f3.2, 

311b  ,  beta]  =■  ,  f  5  . 2  ,  1  1  b  ,  zeta3  *,f5.2/  3x,9h  ’.amma 

3  fO.l,  13b  ,  ai.a  par  *  ,  f  f>  .  3  /  ) 

call  spoctr  (ht,btq,spec,l,aflg,in,gtltl) 
a  f 1 g= .  true. 

w  r  i  t  e  (  6  ,  o  0  3  )  r  0  ,  r  1  ,  r  2  ,  r  3 , 1  p  ,  a  1  f  a 0  ,  a  1  f  a  1  ,  a  1  f  a 2  ,  a  1  f  a 3  ,  be  t  a  0  , 
lbetal ,beta2,beta3,zetaO,zetal , zeta2,zeta3, gamma 
forinat('l  input  new  parameters.  to  enter  a  new  spectrum'/ 

)'  or  aiming  parameter  enter  -1  for  rO '/ 

2 '  i be  old  p  a r  a  me  t  e  r  s  are:' ,  3  i  5 , 4  f  6 . 2  /  4  9  x  ,  4  f  4 . 2  /  4  9  x  ,  4  f  f>  .  2  ,  f  3 . 2  ) 
read  (9,*)  rO,rl tr2,r3,lp,alfaO,alfal ,alfa2,alfaJ,bctaO, 

1  be  t  a  1  ,  be  t  a 2  ,  be  t  a 3  ,  ze  t  at) ,  ze  t  a  1  ,  ze  t  a 2  ,  ze  t  a 3  ,  gamma 

Li  (lp.lt.-l)aflg=.false. 
if  ( lo. It .—1 )lp=— 1 
if  ( l  p  .  g  t  .  1  )  a  f  1  g  =* .  f  a  1  s  e  . 
i  f  (  1  p  .  g  t  .  1  )  1  p  =  l 
i  f  (  ze  t  at) .  1  e  .  0 . 0  )  ze  taO*  1  .0 
if(zetal  .  1  e  .  0  . ')  )  z  e  t  a  1  »  1 . 0 
i  f  (  /.  e  t  a  2  .  1  e  .  0 . 1) )  /.  e  t  a  2  =  1 .0 
i  f  (  ■/.  e  t  a  3 . 1  e  .  0 . 0  )  z  e  l  a  3  *  1 . 0 
i f ( r 0 . e  j . - l ) go  to  30  1 

lp  is  a  switch  with  the  following  meanings: 

lp*0  means  print  only  tiie  summary  table  (brightness  temperature, 
lp»l  means  print  also  the  normalized  contribution  functions. 
li>*  -1  means  only  print  out  the  spectrum. 

lcount’l  con  n  t  1 
go  to  7 


initialize  the  quiet  sun  for  this  a  In par. 

do  4  i  =  1  ,  i  m 

btn(l)“bt(l) 

btq(i)  »  quiet  sun  brig  fitness  temperature  (It), 
go  to  20 


c 
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stop 

e  nd 

/  i  n  c 

s  pek 

/Jut. 

i 

/  i  nc 

s  c  Ii  s  pe  c 

subroutine  s  p  e  c  t  r  (  b  l  ,  h  t  q  ,  s  p  e  c  t  x  l  ,  a  f  1 1; ,  i  in  ,  t  1 1  I  e  ) 
integer  title(15) 

real  *4  1>  t  (  2  1  )  ,  b  t  q  (  2  1  )  ,  spe  c  (  2  1  )  ,  x  1  (  2  l  ) 
r  e  a  1  *  4  q  b  t  (  2  3  )  ,  llmda(23),pspec(23) 
integer  c  li  a  r  /  0/ 

logical  a  f  1  r; 

wrtte(ft  , 9 ')  ) 

9  9  f  o  r  in  a  t  (  '  0  llmda  '  ,  5x,  '  ratio'  ) 

ii  o  1  0  i  =  1  ,  i  m 

i f ( b  t  q ( i  ).eq.0.U)btq(i  )  =  ht(  i  ) 
qbt(  i  )  =  h  t  ( i  ) / b  t  q (  i  )* 100.0 
1  1  n  d  a  (  i  )  =  a  1  o  >;  1  '>  (  x  1  (  i  )  ) 
write  (6,  1  ;  •  0 )  l  1  ni  d  a  (  i  ),ql>t(i  ) 

1)0  format ( '  ' , 2 f  t  0 . 4  ) 

10  continue 
char  =  cliar+l 

if  (char.  L  e  .  11 )  o  to  11 
a  f  1  .  f  al  s  e  . 

w  r  i  1 1>  (  f> ,  1  0  1  ) 

101  format ( '  change  plotter  paper') 
t  1  call  p  1  o  t  s ( 0 , 0 , 1 5 ) 

call  plot(l. 0,1. 0,-3) 

i  f  (  a  f  1 1>  )>o  to  10 

cha  r= 1 

call  s  y  m  h  o 1 ( 1  . 1, 5. 0,0. 14, title, 0.0, 60) 
q  b  t ( im+ 1 ) *  2  0 . 0 
qbt(in+2)=20.0 
do  11  i *  1  , i m 
pspec(t)=spec(  i  )  *  1 00 . 0 
l  1  if(pspec(i).eq.0.0)pspec(i)“100. 
p  s  pe  c (  i  m+  l  )  =  2  0 . 0 
p  s  pe  c  (  i  :n+2  )  =  2  0 . 0 
call  s  c  a  1  e  (  l  1  in  d  a  ,  9  . 0  ,  t  m  ,  1  ) 

call  a  x  i  s  (  0 . 0 , 0 . 0 , 1  0  h  1  o f»  (  1  a .nd a  )  ,  -  1  0 , 9 . 0 , 0 . 0 , 1 1  :nd  a  (  i  n  +  1  )  , 

11  Inula  (im+2)) 

call  a xis(0. 0,0.0, 2 8h brightness  temperature  ratio, 2 3, 5. 0,90.0, 
1 p  <5  p  o  c  ( i  m  + 1 ) , pspec( im+2 ) ) 
call  line(llmda,pspec,im,l,-l,0) 

10  call  line(llmda,qbt,im,l,-l,char) 
call  plot f (0 .0 ,0 .0) 
return 
e  nil 
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Sample  Input  and  Output  at  Computer  Terminal 
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